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INTRODUCTION 



It was the intent of this research to consider the 
potential signal processing use of the charge coupled and 
related devices and how their parameters might help or 
hinder such applications. In order to retain this overall 
perspective, the research was pursued on three levels 
concurrently. Those levels are 

A. Charge coupled device physics . . . the physical 
explanation of the operation of the device in engineering 
terms. We mean to emphasize the parameters the engineer 
is most interested in. A list of such parameters include 
gain, linear range, dynamic range, bandwidth, and maximum 
delay. The determination of the relationship between these 
parameters and the voltages which could be used to clock, 
bias, or input/output the device was paramount. 

B. The CCD as a circuit element. The object of this 
level of the research was to determine how the CCD functioned 
as a sampled-analog delay device. The natural application 
was in delay line filters and both first and second order 
devices were built and tested, relying heavily on the 
characterization performed in A. and also on our growing 
knowledge of filter design, theory and practice. Particu- 
larly emphasized in our research was the recursive comb 
filter since it appeared that in this area we might make 

a significant contribution to technology. 
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C. The application of the circuit element (the recursive 
comb filter) to major electronic systems. System specifica- 
tions were determined to have a profoxand effect on the de- 
vices selected for implementation. It was our goal to seek 
out these specifications and adapt our CCD's (or, in some 
cases charge transfer devices) to the system at hand. The 
most widely known application of the periodic frequency 
transfer f\anction is the radar system moving target indicator 
(MTI) . Physical realities were assessed to evaluate the 
applicability of the CCD filter. 

We believe that this purpose has been, for the most 
part, achieved. Chapter one is the treatment of the param- 
eters of the CCD. Chapter two follows with an exhaustive 
treatment of the filter design problem including the b-plane 
triangle of stability derivation. Finally, chapter three 
deals with the radar system and how requirements dictate, 
or at least, suggest, use of various filter designs. 

Each of the chapters is capable of standing alone. The 
connecting link among them is the use of a saimpled-analog 
device as the primary circuit element. 
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CHAPTER ONE 



THE CHARGE COUPLED DEVICE AND ITS ENGINEERING CHARACTERISTICS 
1 . 0 Introduction . 

Boyle and Smith's paper [1] in 1970 introduced the con- 
cept of causing packets of charge at the surface of a semi- 
conductor to move about on the surface by moving the poten- 
tial minimum (or maximum) . Their brief paper opened the 
door to a new technology which has grown so quickly that 
already three international conferences [73,74,75] on the 
applications of charge coupled devices have been held, each 
opening new frontiers of applications of these devices. The 
reasons for such wide acceptance can be found partially in 
their unique two-dimensional characteristics which permit 
the simple design of imaging devices. This application was 
perceived early, and both theory and experimental results 
were described as early as 1971 [2,3]. 

Another of the attractive characteristics is the sim- 
plicity of fabrication, which rendered the CCD easier to 
manufacture than the older BBD (bucket brigade devices) 

[4]. 

Finally, the charge coupled device is capable of being 
used as a sampled analog circuit device. Unlike most recent 
semiconductor developments which are digital in nature, the 
CCD can and has been used in analog signal processing. In 
particular, the CCD will almost surely replace the older 
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forms of analog signal delay lines (quartz crystals, lumped 
circuit parameters) in many systems since it is 

1. Smaller 

2. Requires no tight mechanical specifications 

3. Potentially cheaper 

4. Less sensitive to temperature variations 

5. Capable of providing longer delay (=1-5 seconds) 

The CCD was first used as the delay element in a trans- 
versal filter in 1972 [5] . Recursive filter design using 
CCD's was not reported until 1974 [6] and by our own research 
at NFS in 1975 [7] . 

Section 1,1 The Charge Coupled Device 

The CCD is an array of MOS (metal-oxide-semiconductor) 
capacitors, each with a separate electrode. These electrodes 
are pulsed in order to move a packet of charge about on the 
substrate. This packet bears a relationship to the informa- 
tion which was used as the input. In case of digital input, 
the packet is present, or not. In the more interesting 
case (for our use) of analog input, it is hoped that the 
size of the packet is related, preferably linearly, to the 
input . 

First we wish to look at the structure of a typical CCD, 
shown in Figure 1.1. Included are the metal electrodes, the 
dielectric oxide, and the n-doped silicon. This implies a 
p-channel device. The (J)'s are the clock pulses. Notice 
that in the configuration in Figure 1.1a, the oxide thicknesses 
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Figure 1.1 CCD structure 
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Figure 1.2 Pulse trains necessary for proper operation of the CCD. 
Figure 1.2a shows THREE PHASE operation. Figure 1.2b shows 
TWO PHASE operation. 
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under each electrode are the same, and that there are three 
separate clock pulse trains (shown in Figure 1.2a) . On the 
other hand, the configuration in Figure 1.1b shows differing 
oxide thicknesses tinder electrodes and only two clock pulse 
trains (shown in Figure 1.2b). For obvious reasons, the 
former is called THREE PHASE operation, the latter TWO 
PHASE operation. The TWO PHASE operation, though more com- 
plicated in fabrication, is widely used because of the sim- 
plicity of clock pulse trains which can drive it. Figure 
1.2 shows each pulse train. In the THREE PHASE operation, 
three distinct levels of voltage are observed. If three 
adjacent electrodes are considered, it becomes clear that 
transfer is effected at the times indicated by a dotted 
line. At the first dotted line, transfer of the charge 
stored under electrode no. 1 to electrode no. 2 is effected. 
At dotted line no. 2, transfer is made from 2 to 3, and so 
forth. This requires the three voltage levels to insure 
that transfer is always made in the correct direction. As 
a matter of fact, this points up an advantage of three 
phase operation in that it can be used to transfer charge 
in any direction and depends only on the magnitude of 
voltages used. As we will see, in TWO PHASE operation, the 
direction of transfer is determined by the structure of 
the oxide thicknesses. 

In the TWO PHASE operation, we see the transfer of charge 
being made when one pulse goes high and the other low. The 
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differences in the oxide thickness enable us to use a two 
level clock pulse because the narrower oxide provides for 
more of the voltage to be dropped IN THE SEMICONDUCTOR 
instead of across the oxide. Thus the surface potential 
(silicon surface) under the two insulator thicknesses are 
different. The surface potential profile is given by 
Figure 1.3. Although this structure does not easily lend 
itself to backward transfer, it is widely used because of 
the simplicity of the clock pulse train. 

We have stated that the CCD is an array of MOS capaci- 
tors. Some of the fundamentals of MOS capacitor character- 
istics are defined here. First of all, if the voltage 
across the MOS device is such to attract majority carriers 
to the oxide-silicon surface, then we say that we are 
operating the device in the accumulation mode . Figure 1 . 4 
shows, at the top, the accumulation mode. 

If, instead, we are operating the device in the reverse 
biased mode, so that the majority carriers in the silicon 
are repelled from the oxide-silicon interface, then we say 
we are operating in either the DEPLETION or the INVERSION 
mode. If the charge (positive) gathered at the interface 
is made up of positive ions, which are relatively immobile, 
then the DEPLETION mode is what we have. If holes appear 
as the positive charged particles, then we have achieved 
INVERSION ... that is, the carriers in the silicon are no 
longer electrons. Instead, the silicon has been inverted 
and the carriers are holes. 
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Figure 1.3 Surface potential, (p- , profile 
of a TWO PHASE CCD. 



ACCUMULATION: 



DEPLETION: 



INVERSION: 




Figure 1.4 MOS capacitor fundamentals 
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Figure 1.5 shows the charge profiles and electric fields 
in the three cases. Notice that the electric field extends 
much further into the silicon for DEPLETION than for accumu- 
lation. This is because the charged region in the silicon 
is much wider in depletion than in acciamulation. As inver- 
sion is reached/ the electric field may one again extend 
only a brief distance into the silicon. Thus we note that 
the electric field is very dependent on the mode of capaci- 
tance we are utilizing. 

Capacitance is defined by 



r - M = 

^ AV d 



( 1 - 1 ) 



where the latter equality is valid for parallel plate capaci- 
tors. Figure 1.6 shows the incremental charge accimulating 
at the interface due to an incremental voltage change. Notice 
that the distance between the charged regions is the same 
after the application of AV as it was before the application. 
This implies that the "d" in the above equation (1-1) 
remains constant. Thus the capacitance is constant. Now 
consider the case in depletion where the AQ forming because 
of the application of AV is made up mostly of positive ions. 
Here there are two different "d"'s no matter what convenient 
method is used for measuring them. Thus we must conclude 
that the capacitance is NOT constant, but that it depends 
on the applied voltage as well as the accumulated charge at 
the interface. Figure 1.7 indicates the variation of the 
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capacitance. 




Figure 1.6a 









1 ^ 1 -^ 

d. 



Figure 1.6b 



Figure 1.6 Incremental charge accumulating at the interfaces of a MOS 
capacitor due to an applied AV. Figure 1.6a shows no change in d, while 
figure 1.6b shows a change in d. Thus, in the capacitor in 1.6b, capacitance 
is a function of applied voltage. 




Figure 1.7 Variation of capacitance in a MOS capacitor. 

For our p-channel device, ^^pppshold ^ negative voltage. 
The ordinate axis appears at an arbitrary voltage. 
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capacitance. Notice that the capacitance is greatest when 
operating in the accumulation mode, and is least when driv- 
ing the capacitor into what is called DEEP DEPLETION. Deep 
depletion occurs when a negative voltage exceeding the 
threshold voltage is impressed on the electrode. For low 
frequency operation, inversion, caused by the thermal gen- 
eration of electrons and holes in the depletion region (dark 
current) is used. But if the MOS capacitor is pulsed very 
fast, the dark current cannot fill the potential "well” 
and thus the capacitance seen is that of deep depletion. 

This is the region in which the CCD operates. Thus we can 
see that the CCD, when used for storage of data, is a very 
volatile memory . . . the absence of a clock pulse train 
renders the device quite useless. 

Figure 1.8 shows the basic model of the CCD. Notice 
that it is two capacitors in series. The following are 
equations used to compute the actual values of this capaci- 
tance. The depletion capacitance is a function of the di- 
electric constants of both the silicon and the oxide, the 
voltage impressed on the electrode, and the doping level of 
the silicon, the oxide capacitance, and the thickness of the 
oxide . 



C 



C 



acc 



(1-2) 
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c = c 

oxide accumulation 




Figure 1.8 Basic model of a CCD: Two capacitors 

in series. 
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where 



SiO- 

q 

Nd 

^si 

th 



V, 



FB 



= relative permittivity of Silicon Dioxide 
= permittivity of free space 

= gate bias voltage 

= charge on electron 

= doping density of silicon 

= relative permittivity of silicon 

= oxide thickness 

= flat band voltage 

= surface potential of silicon 



Equation (1-3) indicates that the equations can be broken 
down even further in that the really important parameter 
is not the voltage on the electrode, but the surface poten 
tial of the silicon. Equation (1-3) gives the relation 
between the gate voltage and the surface potential. These 
equations appear in many references [1] [4]. 

Some useful definitions appear below. 
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(1-4) 



Q 



St 



T 



St I 



dark 




(1-5) 



n = AQ/Q; 



e 



1 - n 



( 1 - 6 ) 



storage time, is defined as the time it takes to com- 

pletely fill the potential well solely by dark current . . . 
that is, by the thermal generation of holes and electrons. 
This number is on the order of seconds which provides us 
with a good bandwidth of frequences which will permit proper 
operation. The full well charge is defined as 



capacitance at deep depletion. 

The second parameter is transfer time. L is the length 
of the electrode and D is the diffusion constant of the 
material. Actually, transfer is complicated and one time 
constant is not enough to characterize completely the 
transfer operation. The first CCD paper ever published, 

Boyle and Smith [1] , indicated that the transfer time con- 
stant was equal to the length of the electrode squared 
divided by four times the diffusion constant. Subsequent 
researchers have found that transfer is actually a combination 




(1-7) 



where C. is the capacitance at inversion and C,, is the 
inv ^ dd 



inv 
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of three different transfer mechanisms and each has a 
different time constant. Also, their time constants change 
with time and distance, so the transfer time is quite com- 
plicated. We shall look at some of the basics of this 
subject later. 

The third parameter is transfer efficiency, or analo- 
gously, transfer inefficiency. Inefficiency is defined as 
the fraction of charge left behind during each transfer. 
Since devices are made which get an efficiency of 99.99%, 
it is a good approximation to say 

Q„ = Q^d - ne) (1-8) 

n o 

where is the magnitude of charge in the zeroth transfer 
gate and is the magnitude of charge in the nth gate. 

If this inefficiency is larger than normal, we may have 
encountered severe noise problems. The main noise problem 
is in the interface states of the silicon-oxide boundary. 
These states trap holes and hold them long enough so that 
they are emitted in a subsequent charge packet. Other 
sources of noise are listed below; 

A. "Edge Effects” 

B. Bulk Traps 

C. Dark Current 

D. Thermal Noise in the input/output circuits 
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The presence of the interface states have led to a method 
of alleviating this source of noise (to first order) by 
simply insuring that these trapping states are filled at 
all times. We do this by what is called the FAT ZERO. In 
order to accomplish the filling of the states at all times, 
we must use, for n-silicon (a p-channel device) , a voltage 
for the HIGH phase of the clock which is near the threshold 
voltage. As we can see, the introduction of charge for all 
time in a potential well means that the "well" is never 
totally empty. Nevertheless, we adjust our thinking to the 
idea that zero is actually some finite size of charge packet. 
Since the trapping states are always filled, the signal is 
passed without distortion (to first order) . 

Section 1.2 Input/Output Considerations 

We will now look at some of the input schemes of the 
CCD. We wish to operate the CCD in the deep depletion 
region, but introduce minority carriers (holes) to the 
electrodes in accordance with some desired signal. The 
first method of input is called DIODE CUTOFF and the physical 
arrangement is shown in Figure 1.9 The signal voltage is 
added to the bias voltage on the input diode . The diode is 
forward biased by this combination of bias and signal. The 
first input gate is pulsed, thus sampling the signal voltage. 
When the gate voltage is high, no charge is injected into 
the n-silicon region. But when the gate is driven low (the 
same level each time) , charge is injected which is related 
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Figure 1.9 DIODE CUTOFF method of input 




Figure 1.10 "Scuppering", "Sloshing" or Potential 
Equilibration mode of input to CCD. V is the signal 
Vpi acts as a constant parametric voltige, and is 
the diode pulse. 
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to the signal voltage in a squared manner. That is, the 
size of the charge packet injected is proportional to the 
square of the signal voltage. Thus we have two disadvantages 
to the diode cutoff method of input; 

1. The pulsed gate voltage must always reach the same 
voltage level. 

2 . We begin our operation in a nonlinear manner . . . 
that is, the very first step toward using the CCD, input, 
is nonlinear. 

For these reasons, we believe that the diode cutoff method 
is not an attractive input mode. 

The BARRIER mode is much like the diode cutoff mode. 

Here the signal "rides" atop a pulse which is applied to the 
diode. The gate voltage is held low. At the arrival of a 
diode pulse, a certain amount of charge is injected into 
the first potential well. It is subsequently transferred 
into the next transfer gate and so on. However, here too, 
we must have precise pulse requirements, and even then, the 
injected charge is theoretically proportional to the square 
of the signal voltage. This nonlinearity makes this mode 
unattractive as well. 

The third mode of input is called POTENTIAL EQUILIBRATION 
or "SLOSHING", or "SCUPPERING". Here the diode pulse is 
used simply as a source of holes to fill the well/wells 
while the signal voltage is applied to one of the input gates 
At the end of the diode pulse, the excess charge under the 
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input gates is drained back across the reverse biased 
p-n junction at the input diode. The residue charge is 
linearly related to the signal voltage at the time of the 
diode pulse. This is best seen by considering Figures 1.10 
1.11 and 1.12. Four instants in time are shown as indicated 
in Figure 1.11. The first time, t^, occurs prior to the 
arrival of the diode pulse. The surface potential profile 
is drawn under the electrodes in Figure 1.12. Time t^^ is 
during the diode pulse. Holes have been injected by the 
forward biased diode into the wells under both and G 2 . 
Although the holes accximulate at the surface, it is possible 
(and very illuminating) to draw them as if they fill UP the 
potential well. After all, when a hole arrives in a region, 
the potential of that region is increased marginally. The 
third time, t 2 f occurs after the end of the diode pulse. 

Since the p-n junction is reverse biased, and since the 
surface potential of the silicon near the junction is very 
negative, the charge drains across the p-n junction leaving 
only the charge which cannot "climb" the potential barrier 
posed by the constant G^ voltage. We have, thus, isolated 
a packet of charge under G 2 which is linearly related to 
the input signal voltage. Time t^ shows the arrangement 
after the low value of arrives. The charge is transferred 
from G 2 to the lower level of the first transfer gate. For 
clarity, the steps of "scuppering" are listed here in order: 

A. Impress signal voltage onto G^^ or G 2 

B. Bias Gj^ and G 2 in accordance with operating curves 
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the input diode pulse 
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Figure 1.11 Time instants relative to the ON/OFF 
cycling of the diode input pulse. 





Figure 1.12a. t=t . Both gates Figure 1.12b. t=t,. Injection of 

1 and 2 are empty. ° holes fills gates 1 and 2. 




Figure 1.12c. t=tp. Excess holes 
drain back across^p-n junction 




Figure 1.1 2d. t=t,. Holes transferred 
to first transfer'^gate when (j), aoes 
1 (>/ ' ' 



Figure 1.12. Sequential action of the CCD in the scuppering input mode. 
Shown :. 3 the stored charge, the input gates, and the surface potential. 
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C. Pulse the diode "ON" briefly 

D. Inject holes into the potential wells under 

and G 2 

E. At end of diode pulse, the holes will move back 
across the p-n diode junction leaving a portion of the 
injected charge. 

F. The magnitude of the residual charge packet is 
LINEARLY related to the magnitude of the signal voltage. 

Figure 1.13 shows the output stage. Notice that there 
exists an output gate designated by V^. This is generally 
lower than any of the clock voltage levels which facilitates 
movement of charge from the last transfer gate, via the 
output gate, and then across the reverse biased junction of 
the output diode. This charge (holes) collects on the 
capacitance at the gate of the output FET, thus turning the 
FET off ... or at least, reducing the current through it. 

At this time, the REFRESH FET is OFF, so that the charge 
simply stays there. When the refresh pulse goes very low, 
this charge is quickly drained off through the refresh FET 
to After draining off, the gate capacitance of the 

output FET is empty of charge and now jusi; awaits the next 
output pulse. The timing of these voltage levels is shown 
in Figure 1.14. 

Figures 1.15 through 1.18 show the input-output charac- 
teristics achieved in our laboratory with a TRW, eight delay 
stage, four input device, by the scuppering input method. 

The segment of the curve on each figure toward the left is 
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Figure 1.14 Timing diagram of the output stage 
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Figure 1.16. CCD operating curves for input through Vg^. Dotted curves show Vg^ as constant 
parametric value, x curves are for Vp^» and A curves are for Vp. 
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the so-called BARRIER mode, which appears to be linear 
despite theory to the effect that it should vary with the 
square of the input voltage. The positive sloped part of 
the characteristic is a noninverting region and is achieved 
by scuppering. The negative sloped part is an inverting 
mode, but it, too, is achieved by scuppering. The explana- 
tion for the shapes of these curves is taken up in the next 
section. 

One significant aspect of CCD research is the question 
of what makes the charge move , . . that is , does it move 
from electric fields or by diffusion. As indicated earlier, 
the first papers assumed that the movement was by diffusion. 
This has subsequently been shown to be incomplete. As a 
matter of fact, there are three phenomena which cause 
movement: 



where is the current density due to diffusion, D is 
the diffusion constant for the material, p is the density 
of holes, and x is the linear distance. 



A. Diffusion 




(1-9) 



B. Drift under the influence of the Fringing Field 



Je^ = q Up P Ef 



( 1 - 10 ) 
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where J_ is the current density due to this drift, u is 

Ef p 

the mobility of holes, q is the charge on one hole, and 
is the fringing electric field. The fringing field is 
the field caused by different surface potentials at different 
regions in the silicon. 



where J_ is the current density due to this drift. 

^s 

The total current is the sum of these currents. Gamma, y, 
in the last equation is a physical constant depending on 
the oxide capacitance (and therefore the dielectric con- 
stant, the thickness of the oxide, and the doping level), 
various voltage levels (gate, flat band, threshold) , and 
other constants. It is interesting to note that all three 
have positive direction at the start of transfer. 

The equation for current density is 



Equation (1-13) is the continuity of charge equation. 



C. Drift under the influence of the self-induced 
field. 




( 1 - 11 ) 



j = qtHpEfP - (YlipP + D)ggl 



( 1 - 12 ) 




9p _ _ 1 aj 

9t ” q 9x 



(1-13) 
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By use of this equation and the expression for the current 
density, we can find a differential equation which relates 
the time rate of change of hole density to the space rate of 
change of hole density and the hole density. This is a 
parabolic differential equation and is not easily solved 
analytically. 

Nevertheless, we can conjecture several facts: First 

of all, at time zero, the dp/dx term is theoretically 
infinite, since in all studies, the initial charge density 
has been assumed to be constant (over x) and zero under the 
gate to which the transfer is about to be made. We conclude 
that the self-induced field is very large, relatively 
speaking. Thus it dominates the current initially. At 
subsequent time intervals, the fringing field, E^, becomes 
the dominant force causing hole movement. As t becomes very 
large, dp/dx actually changes sign and holes diffuse in the 
reverse direction, and tend to cancel the drift current 

caused by E^. When J = 0 or 

/ 

Up E^ p = (Y yp P + D) g (1-14) 

transfer is complete. 

Section 1.3 Experimental Results 

In Figures 1.15 through 1.18, we saw the experimental 
curves relating output to input for our CCD. Given in 
Figure 1.19 are two curves, one indicating the surface 
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Figure 1.19. Gate voltage vs. Surface Potential for the metal and Poly-Si 
gates for the TRW CCD utilized In this research. 



potential of the silicon sxabstrate under the relatively 

O 

thick insulator gate (3000 A) , and the other showing the 
surface potential of the substrate under the thin insulator 

O 

gate (1000 A) . As expected from equation (1-3 ) , the surface 
potential as a function of applied gate voltage is much 
greater for the thin gate than for the thick gate. 

It is possible to explain most of the phenomena shown 
in Figure 1.15 through 1.18 in terms of the shape of the 
curves in Figure 1.19. Such an analysis is important since 
it provides an engineering insight into the physical char- 
acteristics and limitations of the CCD. The analysis is 
perfectly general, since it could be made on any CCD as 
long as the surface potential curves were available. 

The CCD consisted of the following physical devices; 

1. An input diode 

2. Four input gates (of which the second was 
three times as long as the other three) 

3. Eight pairs of transfer gates 

4 . One output gate 

5. One output diode 

6. FET's for refresh and output follower 

The important parameters of this device are as follows: 



O 



1. 


Oxide thickness under 


metal gate . 


. . 3000 A 


2. 


Oxide thickness under 


poly Si gate 


. . 1000 A 


3. 


Permittivity of Silicon Dioxide 


. . 3.9 (relative) 


4. 


Doping level density 


• • • 


1.5 X 10^^/cm^ 


5. 


Metal gate threshold 


• 


. . - 11 V 
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6. Poly Si gate threshold ... - 4 v 

7. Size of input and transfer gates . 1 mil x 1 mil 
(except for input gate no. 2 which is 3 mil x 

1 mil) 

Preliminary Observations (Figures 1.15 through 1.18) 

1. The curves show that there are at least three dis- 
tinct modes of linear operation possible in the CCD, though 
not all present in any given input scheme. For instance, 
if the gate voltage of the gate closest to the input diode 
is varied while another gate voltage is held constant, then 
there is possible only two operating modes. If, on the other 
hand, the gate closest to the input diode is kept at a con- 
stant voltage and some other gate is used as an input for 

a variable signal, then three modes are possible. 

2. The mode with the greatest gain and least input 
linear range of operation is the so-called "barrier" mode. 

3. The modes with the least gain (but still greater than 
one) and widest range are the modes which are usually 
referred to as "scuppering" modes. 

4. The "barrier" mode is an inverting mode. Scuppering 
can produce either inverting or noninverting modes. 

5. There is a local maximum (since the scales are shown 
negative, it is actually a minimum) between the two scuppering 
linear operating regions . 

6 . The barrier mode occurs approximately at the same 
surface potential range for all curves . 
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Explanation of Experimental Results 

Consider Figure 1.15: The gate voltage on is the 

abscissa variable while the dotted curves indicate a con- 
stant parametric voltage on G 2 with G^ and G^ held low. 

The X curves indicate a constant parametric voltage on G^ 
with G 2 and G^ held low. The A curves indicate a constant 
parametric voltage on G^ with G 2 and G^ held low. 

Note that the barrier mode appears to occur for the 
same voltages on G^^ for all input schemes. Since it is 
extremely important to know the actual surface potentials 
under the gates as well as the gate voltages, the Figure 
1.19 is especially important. If we initially set 





^G1 = 


- 15. 


. 0 V 


^G3 = 


- 20 


V 




^G2 = 


- 30 


V 


^G4 = 


- 30 


V 


and 
















^1 = 


- 10 


and - 24 v, 




- 10 


and 


and the 


input 


diode 


pulse : 


12 to - 


2.5 


V 



then the Figure 1.15 shows that we are in empty well 
operation . . . that the voltage on Gj^ is not low enough to 
produce the barrier mode much less the scuppering mode. 
Since we can be sure that the amount of charge transferred 
to the transfer gate is negligible, the surface potential 
diagram shows the situation: (Figure 1.20). The actual 
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G1 G2 G3 G4 T1 T2 




Figure 1.21. Surface potential and stored charge under 
gates 
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values for surface potential in Figure 1.20 and succeeding 
figures have been taken from the curves on Figure 1.19. 



The surface potential for the input diode is taken on 
the n side of the p-n junction. A forward biased p-n 
junction can be expected to have a voltage drop across it 
of about 0.7 V. Thus the -1.8 voltage was obtained by 



When the input diode pulse does not succeed in forward 
biasing the diode (as in this case ) , we do not get a flood 
of holes into the input gate region. However, the jvinction 
IS made less reverse biased, so we do get a few more holes 
injected (by diffusion) than before. Once a hole has crossed 
the junction, it can either re-cross into the p-region, or 
diffuse into the gate regions. We remember that diffusion 
current is established whenever we have a difference in the 
concentrations of carriers in two (or more) different regions. 
Thus, we would expect, for small differences in surface 
potentials of the input diode and to obtain some hole 
injection and subsequent transfer. In analytical terms, the 
current due to the holes is given by 



2.5 - (- 0.7) 



1.8 



(1-15) 




(1-16) 
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Thus as long as 



dp ^ 500 p E 
dx 13 



(1-17) 



we will get hole injection into the gate region. 

When holes reach the region below gate 2, they remain 
there raising the surface potential. They cannot proceed 
further because of the great difference (-20 and -3.5 v) in 
surface potentials between and G^/ and they cannot return 
to for the same reason. They are static charges at this 
point. However, if more and more holes are injected (as 
would be the case if were lowered slightly) , there will 
arise the situation of holes "filling” the potential well of 
G 2 * What this actually means is that the holes attracted 
to the G 2 will finally succeed in raising the surface poten- 
tial of G 2 to about -3.5 v. At this point, some of the 
holes will transfer into the G^ region. Diffusion will 
bring them to G^ where they will remain. Of course, the 
surface potential under G^ will rise at this point. If 
there are enough holes injected during the on time of the 
input pulse, the surface potentials of G 2 and G^ will both 
be -3.5 V and will have charge stored. The surface poten- 
tial of G^ will also be -3.5 v, but will have negligible 
charge stored. This condition is shown in Figure 1.21. 

If still more holes are injected, they will gather under 
gates 2,3, and 4 thus raising the surface potential still 
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further. Finally, at the point of equilibrixom, the surface 
potential vmder all the gates will be -1.8 v. 

When the holes are transferred to the first transfer 
(or delay) gate, goes "low" and goes high. As is 
shown in Figures 1.20 and 1.21, the empty well surface 
potentials are 



high low 

For the metal gate; - 1 v - 4.5 v 

For the poly Si gate; -4.5v -13.5v 

where the fact that the high level of the surface potential 
under the poly Si gate equals the lower level under the 
metal gate is a coincidence. 

The transfer of charge follows precisely the same pro- 
cedures as the transfer from diode to gate. The holes move, 
by diffusion and also under the influence of the fringing 
field and self-induced field, into the transfer gate region. 
The electric field is strong until the holes reach the area 
under the poly-Si. But the deposit of charge in this region 
raises the surface potential. Charge continues to be trans- 
ferred until the surface potential under the poly Si is 
about -4.5 V, the same as the surface potential under the 
metal gate. Thereafter, both regions fill with charge, 
raising the surface potential even higher. The final value 
of the surface potential depends, of course, on the amount 
of charge available for transfer. Thus, if there is a small 
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amount of charge available, then perhaps no charge at all 
will be stored \ander the metal part of the gate. On the 
other hand, if charge is plentiful, then both regions could 
have a final surface potential of -1.8 v. 

Following the first transfer to the surface potential 

diagram might look like Figure 1.22. 

In Figure 1.22 we notice that there still is charge 
stored under even after transfer. This is because the 
surface potentials between and the metal portion of 
are now equal. Thus if holes diffuse into T^^, an opposing 
electric field tends to force them back to G^. A similar 
situation is set up londer G 2 . The charges, permanently 
stored under G 2 cannot be transferred to G^ since the trans- 
fer would set up an opposing electric field which would 
tend to drive the charge back to G 2 . 

Thus far, the mode we have been dealing with has been 
the barrier mode. This does seem to be a descriptive word 
to describe the mode, since the junction is only slightly 
forward biased thus creating a potential barrier to hole 
injection. That this barrier can be overcome is important 
to the operation of this mode. 

When the gate voltage on Gj^ is set lower we immediately 
get full well operation. With the other voltages set as 
before, let us assume that = -17.0 v. Prior to any 
charge transfer, the surface potential diagram is shown in 
Figure 1.23. 
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Figure 1.22 Surface potential and charge distribution 
under gates. 
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Now when the pulse arrives, the junction will be 
strongly forward biased and holes will be injected easily. 

They will first accvunulate under 02^ then , and finally 
bring the entire region to a coitimon surface potential of 
-1.8 V. 

When the input pulse falls, the surface potential at 
the input diode (n-side) will fall to -12.0 v. The charge 
which is stored under will drain through the reverse 

biased diode, leaving G^ empty, which means that the surface 
potential \inder G^^ will return to -2.3 v. Figure 1.24 shows 
the situation at the instant that the input pulse falls. 

In Figure 1.24 it is easy to see that the charge stored 
under G^ will return across the p-n j\inction. It is also 
apparent that this movement will create an electric field 
which attracts the holes \inder gates 2, 3, and 4 to move 
toward the input diode. When the surface potential under 
gates 2, 3, and 4 falls to -2.3 v, the transfer of charge 
across the junction will stop because of the barrier of Gj|^. 
This barrier is set up because at some point in the scuppering 
transfer, further charge movement to G^^ from G 2 will set up 
an opposing electric field ... that is, G^^ will have a- 
higher surface potential than G 2 ^ thus retarding the hole 
flow. After scuppering has been completed, the surface 
potential imder all gates will be -2.3 v, and charge will 
be stored under gates 2, 3, and 4, but NOT under G^. 

The amount of charge transferred to T^ after goes 
low is a fraction of the charge tander gates 2, 3, and 4. 
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gates 
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If the gate voltage of is reduced to lower values, 

LESS charge remains under gates 2, 3, and 4 after the input 
pulse falls. Thus less charge is transferred to when 
(f)j^ goes low. 

Now we set = -20 v. As is shown in Figure 1.15, 

this is near the center of the linear scuppering region. 

Now assume that the input pulse has arrived, fallen to -12 
V, and (|)^ is about to go low. We have the surface potential 
diagram shown in Figure 1.25. 

We notice that no charge is stored under G^. This is 
because the surface potential of an empty G^ is precisely 
-3.5 V. ... the same as it is for an empty G^. When 
goes low, only the portion of charge under G^ is available 
for transfer since the surface potential of G^ represents 
a barrier to hole transfer from G 2 « The figure also shows 
a previously unexplained phenomenon . . . that is, the 
amount of charge transferred levels off for decreased values 
of V_, . This is explained by Figure 1.25. Suppose we 
let V„, fall to -21 v. The surface potential under an empty 
G^, in that case, would be -4 v. The sequence of events 
is shown by the next series of Figures (1.26, 1.27, and 
1.28). In Figure 1.28 we see the final arrangement of 
charge prior to the falling of causing transfer. The 
only charge available for transfer is that under G^. Further, 
we see that an additional decrease in will have no 

effect at all on the amount of charge stored under G^ . This 
amount of charge MUST be enough to cause the surface potential 
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G1 G2 G3 G4 T1 T2 
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G1 G2 G3 G4 T1 T2 




Figure 1.28 Surface potentials and stored charge after 
. input pulse falls and excess charge is drained. 




Figure 1.29 Surface potentials for empty wells. Note that 
the surface potential under G4 is higher than that under Gl. 
This fact will be important. 
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under to be -3.5 v. More charge, of course, cannot 
remain under since it will raise the surface potential 
and cause an electric field to influence its movement left- 
ward. The fact that this is the only charge available for 
transfer, does not, of course, mean that we have full well 
operation. On the contrary, we see that the amount stored 
is not enough for full well operation. 

Moreover, we see that no matter how low the gate voltage 
of Gji^ goes, we can never recover empty well operation with 
the current settings of r ^G4 * 

a leveling off region . . . which is confirmed by the curves. 

A similar phenomenon takes place for variations in 
and a constant parametric voltage, V ^2 (gates 3 and 4 held 
low) . 

The situation of holding as a constant parametric 

voltage and varying is also interesting. The figure 

shows that we first get the barrier mode as explained above, 
and then immediately get scuppering. 

Suppose we set 

= - 19 V = - 7.5 V 

and keep gates 2 and 3 low. Figure 1.29 shows the empty 
well surface potentials. Notice that the surface potential 
under G^ is higher than that under G^. This fact will be 
important . 
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Following the usual chain of events, the input pulse 
causes hole injection which raises the surface potential of 
all gates to about -1.8 v. When the input pulse falls, the 
excess charge is drained off leaving the surface potentials 
of gates 1, 2, and 3 at -3.0 v, and the surface potential 
of gate 4 at -2.8 v. In addition, it is to be noted that 
there is charge stored ONLY under gates 2 and 3, as shown 
in Figure 1.30. 

Now when falls, there is NO charge available for 
transfer since no charge is stored under G^, and there is 
a barrier between the charge stored under and the first 
transfer gate. Thus we get empty well operation. This is 
consistent with the curves in Figure 1.15. 

Now if is raised to -16.5 v, the empty well surface 

potential of G^^ is slightly higher than that of G^. Thus 
after scuppering, we have the following charge and surface 
potential diagram in Figure 1.31. 

A small amount of charge is available for transfer 
according to Figure 1.31. When the charge is transferred, 
the resulting surface potential of gates 2, 3, and 4 can 
be NO LOWER than -2.8 v. Thus, only the charge differential 
between the surface potentials of -2 and -2.8 v can be 
transferred to T^^. 

As is increased, more charge is transferred. 

As is decreased, less charge is transferred. These 

conclusions are consistent with the curves on the Figure 
1.15. 
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If we make the settings 



V 



G2 



5v 



V. 



G3 



and V 



G4 



held low 



and vary the curves tell us that there is a gentle 

slope from empty well operation to almost full well opera- 
tion, We find the surface potentials for this case from 
the Figure 1.19 curves. We notice that the surface poten- 
tial under G 2 is about -1 v. Thus there is a weak electric 
field to inhibit hole movement. The diffusion component 
of current, however, is strong since ^ is large. Thus we 
expect to find a small amount of charge transferred and 
this is, of course, what we get. 

There is one more phenomenon, perhaps not so obvious 
from the curves, which should be considered. It was stated 
that the same phenomena occurred when V ^2 held constant 

as when held constant. That is , there is a leveling 

off effect which cannot be changed by variations in 
There is a certain (depending on the constant value selected 
for Vg 2 amount of charge below which it is impossi- 
ble to go by varying alone. However, the amount of 

charge can be reduced (to zero, it turns out) by reducing 
the constant parametric value of V ^2 ^G3* 

Suppose we set 



V, 



G1 



30 V 



V 



G2 



13 V 



with gates 3 and 4 held low 
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The surface potential (empty well) diagram is shown in 
Figure 1.32. 

In the usual manner, the input pulse injects holes into 
the gate region, raising the surface potential to about 
-1.8 V. Then the pulse falls to -12.0 v and the excess 
charge in the gate region drains across the reverse biased 
p-n junction. We are left with the following charge and 
surface potential diagram (Figure 1.33). 

It would seem, at first glance, that the charge stored 
under gates 3 and 4 would be available for transfer. Such 
is not the case. Notice the surface potential under the 
metal part of the transfer gate which has <p 2 impressed on 
it. This surface potential is -4.5 v. This level is the 
same level as attained by when it goes low. Thus when 
(j)j^ does fall, we have the following surface potential and 
charge diagram (Figure 1.34). 

No charge can be transferred from to T^^ because there 
is already a potential barrier preventing hole injection. 
Thus we get empty well operation as predicted and confirmed 
by the curves . Notice that the condition described is a 
function of the surface potential under the metal gate when 
(j> 2 ^ goes low. This, in turn, is, of course, a function of 
the low level of which for our purposes was kept at 

-24 V. For the above case, if this level were reduced, 
the value of V ^2 necessary for empty well operation would 
also be reduced. In other words, we would get a wider input 
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G1 G2 G3 G4 T1 T2 





-20 -20 

Figure 1.35. Empty well surface potential 
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range and lower gain for this input scheme if the low 
level of could be less than -24 v. 

We now turn our attention to Figure 1.17 which shows 
the variation of '^ith constant ^ Immediately we 

notice some differences between this figure and Figure 1.15. 
Most obviously, we see in two cases, three separate linear 
regions. The first is, of course, the aforementioned 
barrier mode. But the latter two seem to have about the 
same gain, except one inverts and the other does not. There 
is also a local maximum (voltage) between the two latter 
modes . 

Let us bypass the barrier mode as having been explained 
adequately above. Let us make the following settings; 



G1 

G3 



-18 V 
-16 V 



(constant parametric value) 
(variable) 



and put gates 2 and 4 low. 

The initial surface potential diagram is shown in 
Figure 1.35. 

In the usual manner, the input pulse injects holes into 
the gate region thus raising the surface potential in the 
gate region to -1.8 v. When the input pulse falls to -12.0 
V, the excess charge drains through the p-n junction leaving 
the surface potential and charge distribution as shown in 
Figure 1.36. 
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When falls, the charge under is available for 
transfer to If were to be reduced, it would tend 

to reduce the charge transferred since more would be drained 
across the input junction. Thus we have a noninverting 
mode which corresponds to that seen on the figure (dotted 
curve labeled = -18 v) . 

However, there is a limit to the reduction of V _ which 

G J 

can be done and still decrease the charge transfer. If we 
reduce > the empty well surface potential of 

Gji^ is the same as for G^. After the falling of the input 
pulse, this new situation is shown in Figure 1.37. 

This represents the local maximum of charge transfer. 

If is further reduced, say to -19 v, then after scupper- 

ing, there will be THE SAME AMOUNT of charge stored under 
G 2 and G^ as is stored in the local maximum condition, but 
there will be an ADDITIONAL charge stored under G^* And 
this charge is available for transfer. Thus, for further 
reductions in we get MORE, not less charge transfer. 

This corresponds to the curves in Figure 1.17 which shows 
this maximum and inverting/noninverting modes on either 
side. 

Figure 1.38 shows the charge distribution and surface 
potential of the CCD where V ^2 been decreased beyond 

the local maximum. 

Comparisons were made between the actual numbers obtained 
on the curves and the corresponding surface potentials as 
predicted by TRW. Correspondences were quite close. 
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Figure 1.38. Surface potential and charge distribution where 
Vq 2 has been decreased beyond the local maximum 
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Of course, the same explanation can be made for the 
appearance of three modes and a local maximum for varying 

constant parametric values of The numbers are 

a bit different, but, otherwise, the explanation is the 
same. 

One conclusion that we reach is that the various constant 
values of, say should not have any effect on the output 

so long as the maximum has not been reached (by V-,,) . That 
is, the curves, = -18 v vs. = ~20v vs. 

and = -22 v vs. should all be the seime until 

is reduced to the extent of reaching the local maximum. 

This is confirmed by the curves in Figure 1.17. 

The reason this should be so can be seen by considering 
the charge available for transfer . . . it is only the charge 
stored under G^. And this charge is not a function of 
(as long as is low enough to insure that the surface 

potential of is below the surface potential of G^) . 

Figure 1.39 shows a typical charge distribution and surface 
potential diagram for this situation. 

Suppose were decreased. It is clear that this would 
affect the amount of charge stored under G 2 ^ but would have 
no effect on the amount of charge stored under G^. Thus, 
the curves should be common, and it is seen, they are. 

Some of the other parcimeters which we wanted to consider 
are the various values of gain, linear range, and dynamic 
range. Table 1.1 shows the values of gain as a function 
of the gate used as the input gate and the gate used as the 
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constant parameter gate voltages 







G1 


G2 


G3 


G4 


I 

N 


G1 


- 


5.70 


5.79 


4.95 












P 

U 

T 


G2 


2.02 

(-4.4) 


— 


2.09 


2.05 




G3 


1.3 

(-4.14) 


1.3 

(-.269) 


- 


1.33 




G4 


-7.2 


-2.0 


-1.73 


- 



Table 1.1. GAIN (scuppering modes) as a function 
of input gate and constant parameter 
gates 



constant parametric voltage gate (both in scuppering modes) 
Inspection of the table shows that for noninverting scupper 
ing, gain changes little as long as the Scime input gate is 
used. The reasons for this can be found in the same type 
of analysis as was undertaken to explain the shape of the 
operating curves. Consider the noninverting scuppering 
mode of operation with gate 1 serving as the input gate. 

At the time of transfer into the first transfer gate, the 
charge available for transfer lies spread throughout the 
input region: in the potential wells of gates 2, 3, and 

4 (but not gate 1) . The amount of charge in this region 
can, of course, be increased by reducing one of the gate 
voltages (2, 3, or 4) . However this modification changes 
the gain very little because very little of the increased 
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charge finds its way into the first transfer gate. As 
Figure 1.40 shows, the charge in the region under gates 2, 

3, and 4 will increase with a reduction in one of the voltages, 
but the charge available to the transfer gate is only the 
charge which is excess to a surface potential of -4.5 v. By 
excess to a surface potential of x volts, we mean that 
portion of the charge, which, if it were absent, would lower 

the surface potential to x volts. Thus, the available 

/ 

charge changes very little, if at all. 

On the other hand, if the input were to be impressed 
onto gate 2, then the charge available for transfer would 
lie under gates 3 and 4 in the noninverting scuppering mode. 
Figure 1.41 shows a typical charge profile for this situation. 
Since there is less charge available for transfer, we might 
expect that the gain would be smaller, since more voltage 
would be required to draw more charge into the region. But, 
once this situation exists, it is apparent that the actual 
transfer does not depend on HOW the charge was introduced. 

Thus the gain should be the same for all three cases of 
constant parametric voltage. And table 1.1 shows that this 
is correct. 

When we come to input by gate 3, we find that the only 
charge available for transfer lies in the region under gate 

4. Since there is still less charge to be transferred than 
in the first two cases, we expect that the gain would con- 
tinue to decrease, and this is borne out by experimental fact. 
Figure 1.42 shows a typical charge profile for this case. 
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Figure 1.40. Surface potential and charge distribution for a 
typical noninverting scuppering mode with the input impressed 
upon gate 1. Charge available for transfer lies under gates 2, 
3, and 4. The dotted line is the empty well surface potential 
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Figure 1.41. Surface potential and charge distribution for a 
typical noninverting scuppering mode with the input impressed 
upon gate 2. Charge available for transfer lies under gates 
3 and 4. The dotted line is the empty well surface potential. 
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typical noninverting scuppering mode with input on gate 3. 



G1 G2 G3 G4 T1 T2 




Figure 1.44 Surface potential and charge distribution for the 
inverting scuppering mode with input on gate 4, Vq 2 ” parameter 
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There is no noninverting scuppering mode for the gate 
4 input. A little reflection can make this clear. There 
is no way (since these are linearly sequential gates in 
space) to introduce MORE charge into the transfer gate as 
the gate voltage on gate 4 is INCREASED. Increasing the 
gate voltage on 4 reduces the amount of charge under gate 
4 and, more specifically, reduces the amount of charge ex- 
cess to a surface potential of -4.5 volts. Thus reducing 
the gate 4 voltage is the only way to obtain more charge 
transfer. This is a case of inverting gain. As for the 
magnitude of the gain, one can see that it is definitely NOT 
constant. . . that it is a strong function of the gate used 

for the constant parametric voltage. The reasons for this 
can be seen by a look at the charge profile situation 
existing for each case. When is used as a parameter, 

as shown in Figure 1.43, the charge available for transfer 
lies under gates 2, 3, and 4. All the charge in this region 
which is excess to the surface potential under gate 4 is 
available for transfer. Moreover, as Figure 1.44 shows, if 
V ^2 is used as the constant parameter, the charge available 
for transfer is limited to the region under gates 3 and 4. 
Similarly, if is used as a parametric voltage, then the 

charge available for transfer lies only under gate 4. Thus, 
one would expect that the magnitude of the gains would be 
decreasing as the parameter voltages moves from gate 1 to 
3. And this is the result we find in table 1.1. 
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The other inverting scuppering inodes follow similar 
rules except for the strange situation of impressing the 
input onto gate 3 while using V ^2 the parametric voltage. 
The magnitude of the gain is quite small (-0.269) and no 
satisfactory reason has been found for it. One conjecture 
is that the marginal charge available for transfer by 
reducing so small as to be almost negligible. But 

it is not clear why this should be so. The charge available 
for transfer can be shown to lie under gates 2 , 3, and 4 
in this situation, not unlike others we have investigated. 
However, in this case, the CCD behaves abnormally. Such 
a small gain is not very useful for design engineers, and 
thus it probably has no practical significance, yet an 
explanation for its "gain mutation" would be enlightening. 

Other parameters which are important are linear and 
dynamic ranges. The input linear range can be seen from 
Figures 1.15 through 1.18. The widest linear range is about 
7.5 volts on Figure 1.17 using the input and as 

the constant parameter voltage. This is surely adequate 
for most applications. 

The dynamic range for this CCD was approximately 60 db. 

We found that the shape of the output pulse (full well 
operation) of the CCD became considerably poorer as frequency 
was increased. The limit for proper operation was around 75 
kHz, although it could be operated at higher frequencies 
with an attendant loss in gain. We consistently used a 
conservative 20 kHz clock frequency. 
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The maximiim delay time can, of course, be found from a 
calculation of storage time. Empirically, we found that we 
obtained full well operation in all modes at frequencies 
lower than about 12 kHz. This is a result of the dark 
current filling the transfer gates during the "low" time 
of the clock pulses. This corresponds to a storage time of 
83 msec and a total maximum delay time of about half a 
second. Naturally, this calculation is approximate since 
dark current causes considerable distortion even if it does 
not totally fill the well. Better CCD's have been made 
whose delay time is measured in seconds . The watchwords 
of these CCD's are low inefficiency, long registers, and 
a small dark current. 

Section 1.4 Conclusions 

It is important to note that in order to undertake this 
analysis, many voltage levels had to be defined including 
the levels on the refresh FET, clock voltages, the input 
diode bias and pulse levels, plus the output resistance 
(load) value, the width of the input pulses, and the width 
and timing of the refresh pulses. Such a large number of 
parameters is, of necessity, a disadvantage when attempting 
to utilize the device for systems applications. Manufac- 
turers of CCD's must simplify the operation of the CCD such 
that several of the voltages can be tied together at some 
mutually appropriate level. It should be possible by design 
of the geometry and the thicknesses of the oxide levels to 
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insure proper operation of the device with, say 

The user ought to be able to connect the unused input gates 

to as well. 

Some simplicity in timing circuits also ought to be 
designed into the chip. The clock pulse train driving each 
transfer gate is a square wave and the optimum positions 
(in time) for arrival of the input diode and refresh pulse 
are determined, completely, by the transfer clock square 
wave. It is possible in today's high semiconductor tech- 
nology to expect CCD makers to provide one pin for the clock 
and have internal circuitry to provide inversion (for the 
opposite clock phase) and to provide timing signals for the 
input and refresh pulses. 

A serial- in, serial-out CCD designed for scuppering, 
would be left with the following pinouts; 

1. 3. 5. 7. output 

2. Ground 4. G 2 6. G^ 8. clock (master) 

If proper design of the geometry was achieved, and if 
the characteristics of the device could be predicted (as 
we think they can) , then this dream might prove to be a 
reality. As it is at present, there are so many voltages 
and timing problems that it is only with difficulty that 
the devices can be made to operate at all. 

This difficulty was especially troublesome when it came 
to designing and correctly biasing a second order filter with 
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two CCD's serving as the sampled analog delay elements. 

Two CCD's had to be found which were sufficiently "close" 
in their operating characteristics that certain bias voltages 
could be used for each device. This was achieved, but due 
to the drift of the levels, a close watch on the output wave- 
form of the CCD had to be maintained at all times, since 
a drift in the wrong direction of any number of bias and 
clock levels caused the filter to malfunction. 

This last experience brings us to the other urgent 
need: the need for reproduceability . Industry consistently 

tries to sell products whose parameters they can maintain 
within a set of tight specifications. Clearly, the CCD's 
we used in this research (some experimental models from 
TRW SYSTEMS) did not meet this criterion. Nevertheless, as 
will be shown in the next chapter, we were able to design 
both first and second order filters and test their performance. 
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CHAPTER TWO 



UTILIZATION OF THE CCD AS DELAY ELEMENTS IN RECURSIVE FILTERS 
2 . 0 Introduction 

The contents of this section serve to extend the study 
from the CCD itself to how the sampled-analog delay element 
may be used as a circuit element. The following subjects 
will be addressed. 

Section 2.1 . The recursive first order comb filter 

Section 2.2 . Experimental results and theoretical 
calculations of representative first order comb filters. 

Section 2.3 . The recursive second order comb filter 

Section 2.4 . Design methods for a second order 
sampled-analog comb filter. 

Section 2.5 . Experimental results and theoretical 
calculations of representative second order comb filters. 

Section 2.6. Conclusions. 
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SECTION 2.1 . THE SAMPLED ANALOG RECURSIVE FIRST ORDER 
COMB FILTER 



2.I.I . The basic characteristics of a sampled-analog recur- 
sive first order comb filter are straightforward. The filter 
has/ at most, one zero and one pole. The transfer function 
is most conveniently written in the z-domain: 



H(z) 



a + a, z 
o 1 

b + b, z 
o 1 



where 



„n _ jnwT 
z = e 

which has the result that for real frequencies of operation, 
the amplitudes of the transfer function can be found by 
inspecting the functional values on the unit circle of the 
z-plane. By its recursive nature, feedback is present. 
Feedback has a significant effect on the time response, a 
characteristic which plays a major role in applications. 
Recursive filters have the widely assumed bad property of 
having a vmit pulse response which is- infinite in duration. 
Naturally, transversal filters (those without feedback) 
have no such disadvantage. 

Since there is only one pole and one zero, the job 
remains, then, to choose the "best" positions on the z-plane 
for the singularities. With only one pole and one zero. 
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such a task is not especially difficult since the coeffi- 
cients used to set the poles and zeros are numerically re- 
lated to the position of the singularities in a simple 
manner. That is, the design engineer in the laboratory can 
change the potentiometer performing the multiplication 
and know, rather intuitively, what is happening to the pole 
on the z-plane. 

The first method of first order recursive filter design 
consists of simply transforming a known transfer function 
in the s-domain to an "equivalent" one in the z-domain. 

Several types of continuous filters are widely known and 
documented. It is a simple matter to choose a first order 
Butterworth filter:,, for instance, and apply a standard z 
or bilinear z transformation to its s-plane transfer function. 
The resulting function in z is the transfer function of the 
sampled-analog filter. 

As an example, consider the low pass transfer function 
in the continuous time domain: 

H(s) = (2-1) 

1 

If we apply the standard z transform to this transfer func- 
tion, the usual procedure is to find the inverse La Place 
transform (h(t)), and sample it at intervals of, say, T 
seconds. This produces 

h* (n) = h(nT) 
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The Z transform of h' (n) is 



Z(h'(n)) = Z(h(nT)) = Z h(nT)z ^ (2-2a) 

n=0 



1 



A 



-S . T 



e 1 



-1 

z 



(2-2b) 



The form shown in Equation (2-2b) makes it easy to pick off 
the proper values of the coefficients for the sampled-analog 
recursive filter; 

a = 1 b = 1 

o o 

-s . T 

aj^ = 0 bj^ = - e ^ 

where A is the overall gain of the filter. 

The important property of the standard z transform is, 
of course, that the transformation from s-plane to z-plane 
maintains an invariant impulse response. By this is meant 
that the sampled impulse response of the continuous time 
filter is equal to the unit pulse response of the sampled- 
analog filter. Such an invariance is implied by Equation 
(2-2a) . This property is important wherever time domain 
response is a critical specification. 

More often, however, the frequency response transfer 
function is the dominant specification and this is where 
the standard z transform is deficient. The sampling process 
of the standard z transform inevitably produces a finite 
amount of aliasing into frequency response of the sampled-analog 
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filter. The reason for this can be seen by observing that 
the frequency content of the continuous filter is infinite 
... that is, the amplitude of the frequency transfer function 
is nonzero for arbitrarily high frequencies (the only zero 
is at infinity) . When such a signal is sampled at some finite 
sampling frequency, f^, the Nyquist sampling theorem predicts 
that the segment of the transfer function greater than fg/2 
will cause "aliasing” or "foldover". This aliasing produces 
^plitude and phase distortion in the output and prevents 
precise recovery of the original signal no matter what signal 
processing techniques are utilized. This same aliasing in 
the sampled-analog filter causes the amplitude of the transfer 
function to deviate from that of the s-plane transfer func- 
tion at all points along the jw axis. Thus, the frequency 
response of the sampled-analog filter is different from the 
continuous time filter. This result is inevitable and exists 
in all filters derived by the standard z transform. Of 
course, the deviation may be small and the application may 
well tolerate such deviation. But the deviation is always 
present. 

When this deviation becomes large, it is possible for 
important parameters to shift outside the region of accepta- 
bility. For instance, the cutoff frequency^ of a sampled-analog 



^Defined as the frequency for which the response of the 
filter is down 3 db from the response at dc (for low pass 
filters) . 
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filter designed with the standard z transform will always 
be higher than the cutoff frequency of the analog filter 
from which it was derived. In extreme cases (extreme, but 
not unusual) the cutoff frequency may cease to exist at all. 
One example of this is shown in Figure 2-1. The frequency 
response of the sampled-analog filter never takes on a 
value of 0.707H(0). T hus, there is no c utoff frequency. 
Clearly, such a filter, designed for low pass use, is 
worthless . 

The rule of thximb by which sampled-analog filters are 
designed by the standard z transform method is that the 
sampling frequency must be quite large with respect to the 
major frequency content of the continuous- time transfer 
function. To illustrate, we consider Figures 2.2a and 2.2b. 
In Figure 2.2a, the sampling frequency is clearly large 
enough to prevent any significant aliasing from occurring 
in the sampled analog filter since it is much higher than 
the range of frequencies where the transfer function is 
large. However, in Figure 2.2b, the sampling frequency lies 
in the range of frequencies where the transfer function is 
significant. Thus, we would expect substantial distortion 
of both amplitude and phase because of excessive aliasing. 

Another method for finding an "equivalent" sampled 
analog filter for a given continuous time filter is to use 
the bilinear z transform. Here the transformation into the 
z-domain is done directly from the s-domain with no inter- 
mediate step of finding the inverse LaPlace transform of the 
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Figure 2.1. Example of a sampled analog filter whose transfer 
function is such that there is no cutoff frequency. . .no frequency 
for which I H(f) I = .707|H(0)|. 






H'(f) 



Fioure 2.2§ 



Figure 2.2. Examples of the selection of sampling frequency. When 
a large (compared to the frequency transfer function domain) sampling 
frequency is chosen, very little aliasing takes place. However, 
if too low a sampling frequency is chosen, aliasing can make the 
filter useless. 
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analog filter transfer function. Any transformation which 
preserves the requirement that real frequencies be trans- 
formed onto the unit circle of the z-plane will work. How- 
ever, it is customary to use the simplest: 



_ _ 2 z - 1 
® T z + 1 



(2-3) 



Such a substitution requires only algebra to derive the 

transfer function such that the coefficients corresponding 

to a , a, , b , and b, can be found, 
o i o 1 

For our low pass example given for the standard z trans- 
form method, we simply make the siibstitution defined by 
Equation (2-3) and clear of fractions to obtain 



H(z) 



1 + z 



-1 



SiT - 2 
1 



(2-4) 



This form permits us to find the values of the coefficients; 



a 



o 



1 



a 



1 



1 



b = 1 
o 

s.T - 2 

b = 

s.T + 2 
1 



The single most important property of the bilinear z 
transform is that the resulting sampled analog filter assumes 
all the functional values that the analog filter does. That 
is, the range of the functions are precisely equal. Now, 
this does not mean that the equivalent values of the transfer 
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functions occur for the same frequencies. As a matter of 
fact, in general, this is not the case. This change in 
range vs. domain is called "frequency warping" in the 
literature, but it is simply because we are dealing with 
two different functions. Design rules for such filters 
dictate that a desired cutoff frequency in the sampled 
analog domain must be found by solving the following equation 

1 -1 •2>=i 

= 2ir 

1 + b^ 

Conversely, if our cutoff frequency is known (a design 
specification) , then we can calculate b^ by the following 
relation: 

-1 ± sin 2irf T 

b, = ^ (2-6) 

cos 2Trf^T 
c 

As a result of the property which forces the ranges of 
the sampled analog and continuous filters to be equal, we 
must have at least one zero in the frequency range 0 < f < f 
since, as always, the transfer function is periodic. This 
zero can be seen to occur at f_/2, since this is the point 
for which z = = -1. Thus the shape of the transfer 

function is distorted from the shape of the analog filter 
in order to require the zero to occur. This is, of course, 
a good example of what causes the "frequency warping" 
phenomenon . 
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In general/ the relation between the continuous frequency 
variable and the sampled analog frequency variable for which 
equal amplitudes of the transfer functions occur is given by 

= tan -j- (2-7) 

where is the radial frequency in continuous time and 
is the radial frequency in the z-domain. 

Nothing has been said about the unit pulse time response 
of the sampled-analog filter. As a matter of fact, the 
unit pulse time response of the sampled-analog filter does 
not, as in the standard z transform method, equal the sampled 
impulse response of the continuous time filter. For those 
applications where invariance of the time response is a 
necessity, clearly, the bilinear z transform cannot be used. 

The importance, however, of the equivalence in ranges 
of the analog and sampled-analog transfer functions cannot 
be overemphasized. It is this property which affords the 
designer the luxury of obtaining an appropriate analog fil- 
ter (Butterworth , etc.), applying the bilinear z transform, 
and being sure that the amplitude values of the analog 
filter will be faithfully reproduced in the sampled analog 
filter. There is never a question of the existence of a 
cutoff frequency. If it exists in the analog filter, it 
must exist in the sampled analog filter (albeit, at a 
different frequency) . If, as some applications require, 
there is a range of frequencies for which the amplitude 
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response is down, say, 30 db from the maximum response in 
the analog filter, there must exist a range of frequencies 
in the sampled analog domain which exhibits the same behavior. 
This fact has enabled us to emphasize design by bilinear z 
transform in much of our experimental work. 

For reference purposes. Table 2-1 has been included. 

The table shows the relation between the low pass and high 
pass analog filters and their relatives in the sampled analog 
domain. Both bilinear and standard z transforms are per- 
formed on these filters with the resulting coefficients as 
shown. We will explore these effects in greater detail in 
Section 2.2. 

TABLE 2-1 



DESIGN METHOD 


FILTER TYPE 


a 

o 


^1 




Standard z 


Low Pass 


1 


0 


O 

3 

1 

0) 

1 










-o)^T 




High Pass 


1 


0 


c 

e 










0) T - 


Bilinear z 


Low Pass 


1 


1 


- 










0) T ■ 




High Pass 


1 


-1 


+ ^ 

0) T • 



97 



2.1.2 Comb Filters 



The basic characteristics and the design rules of first 
order filters are so simple that it is apparent that anytime 
a filter of higher order can be put into the first order 
form, it is a wise thing to do. For instance, if we have a 
sampled- analog filter of, say, eighth order, and if it can 
be shown to be equivalent to some first order filter, then, 
we need only to consider the first order filter, and, at 
the conclusion of our design, make the proper translation 
into eighth order. 

In our research, we have used a serial-in, serial-out 
charge coupled device as our analog delay element. The 
specific chip we have used contained eight delay stages so 
that eight samples of the input were stored in the CCD at 
all times. The filter which we designed utilizing this 
analog delay line consisted of multiplicative coefficients 
at the input and output of the CCD, followed, of course, by 
high input impedance summers. The filter set-up is shown in 
Figure 2.3. Notice that the coefficients, a^ and bg, 
multiply the output voltage by an appropriate value prior 
to svunmation. The transfer function of this filter is 



H(z) 



1 + bgz‘® 



(2-8) 



The important characteristic of this transfer function is 
that the coefficients, a^^,a 2 , •••, and bj^,b 2 ^ .«., 
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input 




Figure 2.3. Experimental set-up for first order 
recursive filter 




An example of the eight to first order translation 
possible with serial-in, serial-out delay elements. 
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are all zero. Clearly, this is an eighth order filter. 

It has eight poles and eight zeros. In order to find the 
position of these poles and zeros on the z-plane, it is 
necessary to solve the equations 

z® = and z® = -bo (2-9) 

^o ® 

Such a solution involves the eighth root of some real 
number (the coefficients, of course, cannot have nonzero 
imaginary parts) . Now, poles lie at the various positions 
of the eight eighth roots of ~bg* In accordance with the 
mathematics of complex variables, each of these roots will 
have equal absolute values or moduli. The angles will, in 
general, be different. The mathematics of complex variables 
also tells us that the angles between the adjacent roots 
will be equal. That is, the poles {and zeros) will be 
EQUALLY SPACED angularly throughout the z plane. Figure 
2-4 shows an example of such a plot. 

Kuo [(H)] has shown that, in the same manner as is 
done in analog system transfer functions, the response of 
the sampled analog filter in the frequency domain can be 
obtained by graphical means. This method consists of picking 
a point on the unit circle and determining the length of 
the line segments drawn to each singular point on the plane. 
The lengths of these line segments are arranged so that the 
following quotient is obtained; 
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( 2 - 10 ) 



jj _ distances from point to zeros 
distances from point to poles 

and H is the absolute value of the frequency response of 
the filter at the frequency of interest (determined by the 
point on the unit circle) . 

Using this method, we can investigate the frequency 
response of our eight pole-eight zero filter. From dc to 
one eighth of the sampling frequency, f /8, the distances 

o 

to the singular points are found and the expression in 
Equation (2-10) is evaluated. In this way a point-by-point 
computation of the transfer function is made. When we begin 
to consider values of f larger than f /8, however, we notice 

o 

that the singular points are in precisely the same relative 
location as they were for appropriate frequencies between 
dc and f„/8. Thus, it is clear that the distances will be 
exactly equal, and further, the expression for H will be 
exactly equal for frequencies f, f+f /8, f + 2f /8, ..., 
f + nfg/8. 

This observation, in turn, is the very nature of periodic- 
ity. Thus, the transfer function values are periodic with 
respect to fg/8. 

To backtrack a little, it is well known that a sampled_ 

analog filter is periodic with respect to f and is symmetri- 

cal cibout f„/2. What we have done is to determine a whole 
s 

class of recursive sampled analog filters which are not only 
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periodic with respect to f , but also periodic with respect 

to fg/8. It turns out that this also implies that our new 

filter is symmetric with respect to f /16. 

s 

The use we are going to make of this discovery was fore- 
shadowed by the first paragraph in this section. This filter, 
although nominally eighth order, can be considered a first 
order recursive sampled analog filter with an adjusted 
sampling frequency of f_/8. Of course, certain adjustments 
must be made in the coefficients as well. For instance, in 
our set-up shown in Figure 2.3, if bg = 0.7, then the poles 
of the transfer function must lie on the ^77 circle on the 
z-plane. However, when the translation to first order is 
made, the pole formed by bg will lie on the real axis (since 
there is only one pole) and its modulus will be 0.7. 

It should be understood, that for clarity, we have used 
in this discussion an eighth order filter. The theory and 
proof are equally good for any order filter as long as the 
condition that the middle coefficients are zero is met. An 
analytical proof is given in Appendix A. 

The improvement in the actual design procedure by using 
this translation is perhaps better appreciated if one vis- 
ualizes the design engineer at his bench attempting to find 
the best coefficient for a given ideal transfer function. 

Of course, there may be other methods (such as with a com- 
puter) which he could use, but one time-honored method is 
actual laboratory experiments. Unless he is unusually 
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bright, he cannot project the changes in the amplitude of 
the transfer function based on a small change of a coeffi- 
cient in an eighth order filter. On the other hand, we 
have seen how simple it is to project changes in the trans- 
fer function of a first order filter by virtue of changes 
in the single coefficient, Thus the engineer is pro- 

vided with an intuitive aid to design . . . without resort to 
the digital computer. 

A more fundamental reason for design in this manner con- 
cerns the construction of comb filters. Comb filters will 
be defined here as those electronic devices which exhibit 
a periodic transfer function, and which are little affected 
by the various nonideal circxjmstances which inevitably ap- 
pear in any physical system. For instance, our theory has 
depended on the ability to sample a waveform instantaneously 
and store a number which corresponds to the value of the 
waveform at that instant. In addition, we have assumed that 
we can restore the original signal by using an ideal low 
pass filter, when, in fact, most of the restoration of 
signals is done with a sample-and-hold device, which is NOT 
ideal. Finally, we have assumed that our delay line is free 
from attenuation and distortion. In a charge coupled device, 
we always have a certain amount of charge left behind during 
each transfer from gate to gate. The charge left behind adds 
to the succeeding packet of charge and causes amplitude 
distortion. These factors, plus other instrumentation 



103 



realities, prevent us from ever producing a truly periodic 
transfer function. The sample-and-hold effect causes a 
sin x/x envelope to be superimposed in a multiplicative way 
on the samp led- analog filter transfer function. The non- 
zero sampling time also causes an additional six x/x envelope 
to be superimposed multiplicatively. The order of these two 
effects is not the same, however. The six x/x caused by the 
sample-and-hold is far more restrictive since the first 
zero of the sin x/x function occurs at f = f^. The first 
zero of the other effect occurs at f = 1/t where t is the 
sampling time width, usually quite small. 

Thus, in designing a comb filter by conventional means, 

the use of sample-and-hold circuitry inherently causes a 

rather rapid fall off of amplitude response. That is to 

say, if there is only one maximum of the transfer function 

between dc and f^, then the second maximum, located between 

f^ and 2f^, is much smaller than the first. The ideal comb 
s s 

filter, of course, has equal maxima throughout the spectrum. 

Now consider the use of the method explained above. 

There are as many maxima between dc and f^ as there are 
delay elements in the delay line, which is, of course, equal 
to the true order of the filter. In our case, this number 
was eight. Thus eight maxima appeared in the region between 
dc and f^. The sample-and-hold circuitry still caused the 
sin x/ X rolloff of amplitude response, but its effect was 
negligible for most applications for the first two maxima. 

The real payoff comes when additional delays are built into 
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the system. One such circuit which we have used is the 
RETICON SAD (serial-analog delay device) . This circuit has 
96 delay elements and the sampled analog filter transfer 
fvinction (using this device) has, in theory and practice, 

96 maxima between dc and f . Clearly, the first twenty or 
so are not affected adversely by the sin x/ x factor. 

Thus, it is possible, using these techniques, to design 
an almost ideal comb filter while, at the same time, offer- 
ing to the design engineer the enormous advantage of being 
able to treat his filter as if it were a first order filter 
with the accruing simplicity of design. Figure 2.5 shows 
the first four "teeth" or periods of a comb filter designed 
by conventional means, corrected for the effect of the sample- 
and-hold circuitry. Figure 2.5b shows four "teeth" or periods 
of a comb filter designed by the aforementioned technique 
using a delay of eight stages. Finally, Figure 2.5c shows 
the same four "teeth" using a 96 stage delay device. The 
comparison between the three is clear. The latter exhibits 
the best performance for periodicity. 

SECTION 2.2 EXPERIMENTAL RESULTS AND THEORETICAL CALCULATIONS 
OF REPRESENTATIVE FIRST ORDER RECURSIVE 
COMB FILTERS 

2.2.1 With the theory of Section 2.1 in mind, the next step 
was to determine whether adequate first order filters could 
be implemented using the charge coupled device as the delay 
element. The CCD selected for the experiment was built by 
TRW SYSTEMS and consisted of an eight bit serial-in. 
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serial-out delay register. The cross section of the device 
is shown in Figure 2.6. This is the same device as was 
used for the investigation of the theory of CCD operation. 

The curves of operation derived in Chapter one were uti- 
lized here to find an appropriate biasing point. To choose 
the optimum point, we desired to have the widest possible 
linear range for both input and output. Such a point was 
found to occur whenever the first two input gates (G^ and 
G 2 ) were biased low (-30 v) , the input impressed upon gate 
#3, while using the voltage on gate #4 as a parameter. The 
actual operating curve is shown in Figure 2.7. The linear 
ranges are easily seen. The steep slope on the left is the 
so-called barrier mode of operation which was discussed in 
Chpater one. The second region is the region of "scuppering”, 
"sloshing", or "potential equilibration". The actual 
quiescent point used was 

= - 18.7 V 

V, = - 12.0 V 
4 

which provided a linear range of approximately 3.7 v. 

Figure 2.3 shows the experimental set-up. For a^ = 1, 
a direct connection was made from the input summer to the 
output summer. For a^^ and b^, potentiometers were connected 
into the circuit after the sample-and-hold and level shifter. 
The full value of these potentiometers was 20 Kfl. The value 
for a given coefficient was determined by considering the 
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Figure 2.6. 



Cross section of CCD. 




Figure 2.7. Operating curves for CCD selected 
for use in filter design. 
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gain of the CCD eind the input impedance of the summers . 

The value of the potentiometers was, then, not obviously 
related to the value of the coefficients in the transfer 
function. Specifically, if the input impedance of the 
summer could be neglected, then, the value of the potentio- 
meter, P, was related to the coefficient, a, by 

P = a/K 



where K is the overall gain of the CCD, sample-and-hold 
device, and level shifter. 

Unfortunately, we found, for first order experiments 
that the operational amplifier summers did not have large 
enough input impedances to justify neglecting them. Thus, 
we needed a correction factor for the stammers. The applica- 
ble circuit diagram is shown in Figure 2.8. The input to 
the circuit (output circuitry of the level shifter) is 
and the desired output is Using the usual virtual 

ground assumption for operational amplifiers, we can calcu- 
late our desired value of the potentiometer which yields 

aV. . We call this value P. The total impedance as seen 
in 

from the output of the level shifter is 



Z.. = 20K - P + 
tot 



PZ . 

in 

P + Z . 

in 



( 2 - 11 ) 



109 




kV 



in 




Figure 2.8. Block diagram of implementation of feedback and 
feed forward coefficients. 




Figure 2.9 . Voltage divider circuit for coefficient 
implementation. Correction factors were calculated which 
included the effects on was very small. 
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This becomes a voltage divider network as shown in Figure 
2.9. Thus our voltage divider transfer fxinction becomes 



aV. 



in 



k 



V. 

in 



P 



Z . 
in 



P 20K + 20K Z. 

in 



( 2 - 12 ) 



In the above equation, all the quantities except P are 
known. Thus, in order to assist us in our experiments, a 
correction curve was generated which related "a" to "P" 
for all appropriate values of a (a < 1) . 

In order to have a coherent set of experiments, it was 
decided to utilize both the standard z and bilinear z 
transforms to generate first order filters. In addition, 
we required that a number of integrator and canceller types 
of filter be investigated. The integrator type filter is 
defined as having narrow pass bands (high Q) with wide stop 
bands where attenuation is high. The canceller is the in- 
verse of the integrator. The canceller has narrow stop- 
bands (high Q where Q for stopbands is defined as the Q 
for a filter whose transfer function is 1 - H(f)) with wide 
passbands. These two types of filters are shown in Figure 
2 . 10 . 

Naturally, it was desired to look at low pass and high 
pass type filters as well. Using these designations, it 
is possible to consider a low pass, integrator type, filter 
designed by the standard z- transform or any other combination 
of these three categories as shown in Table 2.2. 
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Figure 2.10 Frequency characteristics of 
t.vo general types or comb filters - - - 
bands top, canceller type; and bandpass, 
integrator type 
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PASS 


TYPE 


DESIGN METHOD 


LOW 


INTEGRATOR 


STANDARD Z 


HIGH 


CANCELLER 


BILINEAR Z 






OTHER 



Table 2.2: Categories of first order filters 

In order to have general expressions for the coefficients 
(so the coefficient would not have to be derived every time 
we needed a change in its value ) , we made the following 
calculations : 

For the standard z transform, we assume a form for the 
low pass filter to be 

0) 

H(s) = 

s + u 

X 

and find the standard z transform of this H(s). 
is (except for a gain factor) : 



-w T , 

1 - e ^ 

Such a filter is an integrator if the cutoff frequency is 
less than 0.25 times the sampling frequency, and a canceller 
if the cutoff frequency is greater than 0.25 times the 
.sampling frequency. These statements are the logical results 
of our definition of integrator and canceller. When the 



(2-13) 
The result 
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cutoff frequency is exactly 0.25f , the passband is pre- 

s 

cisely as wide as the stopband and neither an integrator nor 

a canceller can be achieved. On the other hand, it is well 

to keep in mind that, by using the standard z transform, 

we may not get any cutoff frequency at all. 

The requirement for the integrator is that the cutoff 

frequency be less than 0.25f (0.25/T). This translates to 

-co^T 

a requirement that b^^ = -e ^ be less than -0.268 (this 
calculation is found in Appendix B) . A quick look at 
exponential tables tells us that, in order to have a low 
pass integrator, we must have f T less than 0.21. This 
last statement tells us something about the relationship 
between the s-plane analog filter and the sampled analog 
filter of the z-plane. The o) is the radial cutoff fre- 
quency of the analog filter. If we desire a narrow band 
low pass filter in the z-plane, we must choose a sampling 
frequency, f , such that f < 0.21f or a sampling frequency 

o ^ S 

considerably higher than the analog cutoff frequency. Sup- 
pose that the analog cutoff frequency was 10 kHz. Then in 
order to design a low pass integrator, we must require a 
sampling frequency of at least 

10 ^ 

's = ^ <2-15) 

Failure to do so may yield a wide band low pass filter, or, 
worse, a quasi all pass filter (one with no cutoff frequency). 
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Naturally, to achieve a low pass canceller, we need to 

have f /f >0.21. We also have a lower limit below which 
X s 

we get no cutoff frequency. This is foiind by determining 
at what sampling frequency the cutoff frequency no longer 
exists. This calculation is found in Appendix D. The 
critical value for b^^ such that if bj^ takes on any larger 
value, no cutoff frequency exists is -0.172. This translates, 
via exponential tables to 

- w^T = - 1.76 (2-16) 

which further reduces to the requirement that 

< 0.28 (2-17) 

Thus , we also have determined what the appropriate 
sampling frequency will be in order to design a low pass 
canceller: 



0.21 < fyf^ < 0.28 (2-18) 

Similar calculations can be made for high pass filters. 
Except for the sign of b^ (it becomes positive since the 
poles must lie on the negative real axis of the z-plane) , 
the values are the same: 
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INTEGRATOR: 



0 < f /f <0.21 
^ s 



CANCELLER: 0.21 <f/f^ < 0.28 

One other possibility exists in using these standard 
z trcinsforin filters. If we simply subtract the transfer 
function (normalized to one) from unity for all frequencies, 
we can achieve the inverse of each filter type with different 
coefficients. For instance, if we subtract the transfer 
f\anction for the low pass filter 



H(Z) = 


tt (2-19) 

1 - h^z 



from one, we get 



H(z) = 


-b.z"^ 

T- (2-20) 

1 - b^z" 



which is a high pass filter. If the original filter was an 



integrator type, then the 


latter filter is a canceller. 


This translation produces 


the following results 


Start with 


TRANSLATION (l-H(z)) Result 


Low Pass integrator 


High pass canceller 


Low pass canceller 


High pass integrator 


High pass integrator 


Low pass canceller 


High pass canceller 


Low pass integrator 
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One remark about these filters is worth making here. 

We observe that for all these filters, only one term appears 
in the numerator, either a^ or a.^z in either case, it 

is clear that no true zero will ever appear in the transfer 
fxinction for real frequencies. Looking at the pole-zero 
plot, this means that a zero will never appear on the unit 
circle. If we are seeking a large passband-to-stopband 
attenuation ratio, it appears that the only means to achieve 
our goal is by seeking the optimum value of bj^. Since b^^ 
is a fxinction of the ratio of the analog cutoff frequency 
to the sampling frequency, this meeins that in searching for 
optimxim values of b^^, we are actually searching for the 
optimum ratio of analog cutoff frequency to sampling fre- 
quency. To be more explicit about this procedure, let us 
consider an example; If we have transformed an analog 

filter whose cutoff frequency was f using a sampling fre- 

-u T 

quency of f^, we have obtained a b^^ of ±e . Now, if we 
observe the frequency response of this filter and decide 
to shift it by changing b^^ slightly, what we have actually 
done (assuming our sampling frequency remains constant) 
is to change the original analog filter from which our 
sampled analog filter was derived. 

Tcible 2.3 shows the generalized formulas for the 
coefficients for a filter derived by the standard z transform. 

Figures 2.11 through 2.14 show the experimental data and 
the theoretical shape of the transfer functions for four 
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TABLE 2-3 



Design Summary of First Order Filters 
Using the Standard z Transform 



FILTER TYPE 



Low Pass Integrator 

fc 

^ < 1/4 
s 



-e 



-0) T 
c 



Low Pass Canceller 



-0) T -0) T 

c c 

-e e 



> 1/4 



High Pass Canceller 



-0) T -0) T 

c c 

-e -e 



^ < 1/4 
s 



High Pass Integrator 



-(ji T 
c 



> 1/4 



X 
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Figure 2.11, Theoretical and Experimental transfer function 
for a first order recursive filter, a = 1 , = 0, b, = 0.3 

Experimental data is shown bv a (+). ° ' * 
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Figure 2.12 Theoretical and Experimental transfer function 
for a first order recursive filter, a =l,a, = 0, b. = 0.7 
Experimental data is shown by a (+). ° ' 
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Figure 2.13 Theoretical and Experimental transfer function 
for a first order recursive filter, a = 1, a, = 0, b = -0.3 
Experimental data is shown by a (+). ° ' 
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Figure 2.14. Theoretical and Experimental transfer function 
for a first - order recursive; filter, a = 1, a = 0, b-j ~ -0.6 
Experimental data is shown by a (+). ® 
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filters. The first filter was designed as a high pass 
integrator. The feedback coefficient is fairly small, 

0.3, which corresponds to a t/t ratio of 0.19. From the 
discussion above, it can be confirmed that this ratio falls 
into the range of a high pass integrator. Nevertheless, 
the ratio is "close” to the limit, 0.21, and therefore, we 
would expect that the filter would have a poor quality 
factor, where Q is defined by 



Q = 



f 



center 



BW 



( 2 - 21 ) 



where BW is the bandwidth of the filter. 

Increasing the feedback coefficient, b^, to 0.7 enables 
us to design and test a filter with considerably higher 
selectivity (a higher Q) . Figure 2.12 shows the results 
of the theory and experiment. 

Figure 2.13 shows our results for small feedback 
(bj^ = -0.3) and a low pass filter. The shape of the fre- 
quency characteristic is the same as for the high pass case 
with the same absolute value of feedback coefficient. 
Therefore, it is, of course, a low pass integrator. 

Finally, Figure 2.14 shows the results for large negative 
feedback. Apparently, during the experiment, the feedback 
coefficient was less than had been expected since the 
experimental and theoretical curves do not correspond so 
closely as the other three cases. 
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We observe that, as expected, none of the four cases 
show a true zero either experimentally or theoretically. 

This is in accordance with our predictions. Agreement 
between theory and practice was good except for the condi- 
tion of large negative feedback. These curves all show a 
linear plot, so that the maximum attenuation ratio between 
peak and null occurs for the higher feedback case and is 
approximately 12 db. The same ratio for the condition of 
low feedback (b^ = 0.3) is less than 6 db. Clearly filters 
should be able to achieve better passband-to-stopband ratios 
than this. A search for better filters leads us to the 
next section. 

2.2.2 The Bilinear Z Transform 

We start with the same analog filter as before and 
apply the bilinear z transform directly to the transfer 
function. 



H (s) 






X 



S + OJ. 



Applying the bilinear z transformation; 



2 z - 1 
T z + 1 



we get 



H(z) 






-1 



{(O^T + 2) (1 + 




03 T - 2 

03xT t 2^ 



( 2 - 22 ) 



(2-23) 



(2-24) 
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Now using the same criterion for integrator and can- 
celler type filters, we must find the applicable ratios of 
f /f_ and appropriate values of a , a., and b, to use. 

One good feature of the bilinear z transform is that it 
produces both nonzero a^ and a^^. In this way, we can be 
sure that a theoretical zero will occur in the transfer 
function. Of course, we will not expect an infinite 
passband-to-stopband response ratio, but the ratio should 
be better than in the standard z transform case. We note 
that the value for a^ and a^ are precisely equal in Equation 
(2-24) . On the z-plane this corresponds to a zero on the 
unit circle at z = -1.0. Since the transfer function will 
be continuous, we know that there will be a cutoff frequency 
. . . there is no uncertainty as there is when using the 
standard z transform. Of course, the cutoff frequency will, 
in general, not be the same cutoff frequency as existed in 
the analog filter since the two functions of frequency are 
different. But a cutoff frequency is certain to exist. 

Requiring that the cutoff frequency of the sampled 
analog filter be less than 0.25f , we find the following 
relations: 

A. For the low pass integrator, the constraints 
are 

-1 < b, < 0 and 0 < t /f. < l/ir 

X X s 
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F. 





j^TP 




B. For the low pass canceller, the constraints 
are 

0 < b, <1 and 1/ir < f /f < » 

where the value of the ratio rises without bound as 

b^ tends toward 1. The physical explanation for this is 
that the cutoff frequency in the analog domain may be as 
large as desired, and still, there will be a cutoff fre- 
quency in the sampled analog domain for ANY sampling fre- 
quency. Therefore, suppose we have a filter with a large 
cutoff frequency, say 100 kHz. Now suppose we wish to apply 
the bilinear z transformation to this transfer function and 
use the resulting filter in a system where the sampling 
frequency is 20 kHz. Now this sort of thing simply could 
not be done using the standard z transform, since, as 
explained in the preceding section, the sampling frequency 
must be considerably larger than the analog cutoff frequency 
in order to achieve proper cutoff in the sampled analog 
domain. However, using the bilinear z transform, we can 
be assured that the cutoff frequency will exist, and it will 
be greater than 0.25fg or 5 kHz and be between 5 kHz and 
10 kHz for a canceller. The power of the bilinear z trans- 
form is emphasized by these computations. 

As a corollary to the existence of the cutoff frequency, 
we can also infer the existence of any other value of the 
transfer function so long as the value existed in the analog 
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filter transfer function. This occurs because the bilinear 
2 transform maps the entire jco axis of the s-plane onto 
one revolution of the unit circle in the z-plane. So, just 
as we always have cutoff, we also always have a point which 
is, say, 30 db down from the maximum, as long as such a 
point existed in the analog filter function. 

In the high pass case, the transfer function derived 
by the bilinear z transform is almost the same as for low 
pass except the cpefficient a^^ is negative. Thus, this 
implies the existence of a zero in the z-plane at z = 1 
(or f = 0, dc) . The coefficient, b^, once again varies from 
-1 to 1 as follows; 

A. For the high pass canceller, the constraints are 

-1 < b, <0 and 0 < f/f < 1/tt 

X ^ s 

B. For the high pass integrator, the constraints are 

0 < b, <1 and <» > t/t > 1/ir 

^ ^ s 

Table 2.4 shows the generalized formulas for the coefficients 
for a filter derived by the bilinear z transformation. 

Figures 2.15 through 2.28 show various filter transfer 
fxinctions for various coefficients. They are tabulated 
here for easy reference: 



Figure no. 


^o 


^1 


®1 


Type 


2.15 


1 


-1 


0.9 


High pass integrator 
(high Q) 



TABLE 2-4 



Design Summary of First Order Filters 
Using the Bilinear z Transform 



FILTER TYPE 


a a, 

oil 


Low Pass Integrator 


1 1 -b^ 


fc 

^ < 1/4 
s 




Low Pass Canceller 


1 1 +b^ 


fc 

^ > 1/4 
s 




High Pass Canceller 


1 -1 -b. 



fc 

^ < 1/4 



High Pass Integrator 


1 -1 +b^ 


fc 

^ > 1/4 
s 




where, in all cases, 


0) T - 2 

u _ c 

°1 w T + 2 

c 
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Figure 2,15, Theoretical and experimental transfer function for 
a first order recursive filter, a = 1, a = -1, b = 0.9 
Experimental data are shown by a f+). ^ ' 
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Figure 2.16. Theoretical and experimental transfer function for 
a first order recursive filter, a = 1, a = -1, b = 0.7 
Experimental data are shown by a f+). ' 
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Figure 2.17. Theoretical and experimental transfer function for 
a first order recursive filter, a = 1, a, = -1, b = 0.3 
Experimental data are shown by a ?+). ' 1 
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Figure 2.18. Theoretical and experimental transfer function for 
a first order transversal filter, a = 1, a = -1, b = 0 
Experimental data are shown by a (+/. 
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Figure 2.19, TTieoretical and Experimental transfer function for 
a recursive first order filter, = 1 , a, = -1 , bi = -0. 3 
Experimental data are shown by a f+) ' 
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Figure 2.20. Theoretical and experimental transfer function for a 
first order recursive filter, a = 1, a = -1, b = -0.7 
Experimental data are shown by 9 (+). 1 ' 
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Figure 2.21. Theoretical and experimental transfer function for a 
first order recursive filter, a =l,a =-1,b =-0.9 

Experimental data are shown by § (+). ' ‘ 
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Figure 2.22. Theoretical and experimental transfer function for 
a first order recursive filter, a = 1, a. = 1, b = -0.9 
Experimental data are shown by a ®+). ' ^ 
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Figure 2.23. Theoretical and experimental transfer function for 
a first order recursive filter, a = 1, a = 1, b = -0.6 
Experimental data are shovm by a ?+). ‘ ' 
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Figure 2.24. Theoretical and experimental transfer function for 
a first order recursive filter, a = 1 , a, = 1 , b, = -0. 3 
Experimental data are sfiown by a '^+). 
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Figure 2.25. Theoretical and experimental transfer function 
for a first order transversal filter, a =1 , a = 1, b =0 
Experimental data are shown by a (+). ° 1 ^ 
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Figure 2.26. Theoretical and experimental transfer function for 
a first order recursive filter, a = 1, a, = 1, b = 0.3 
Experimental data are shown by a ?+) . 
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Figure 2.27. Theoretical and experimental transfer function for 
a first order recursive filter, a = 1, a = 1, b = 0.7 
Experimental data are shown by a f+). ^ ' 
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Figure 2,28. Theoretical and experimental transfer function for 
a first order recursive filter, a = 1, a^ = 1, b = 0.9 
Experimental data are shown by a ?+). 
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Figure no. 


A 

o 


^1 


®1 


Type 


2.16 


1 


-1 


0.7 


High pass integrator 
(medium Q) 


2.17 


1 


-1 


0.3 


High Pass integrator 
(low Q) 


2.18 


1 


-1 


0.0 


High Pass (nonrecursive 
filter) 


2.19 


1 


-1 


-0.3 


High pass canceller 
(low Q) 


2.20 


1 


-1 


-0.7 


High pass canceller 
(medium Q) 


2.21 


1 


-1 


-0.9 


High pass canceller 
(high Q) 


2.22 


1 


1 


-0.9 


Low pass integrator 
(high Q) 


2.23 


1 


1 


-0.6 


Low pass integrator 
(medium Q) 


2.24 


1 


1 


-0.3 


Low pass integrator 
(low Q) 


2.25 


1 


1 


0.0 


Low pass (nonrecursive 
filter) 


2.26 


1 


1 


0.3 


Low pass canceller 
(low Q) 


2.27 


1 


1 


0.7 


Low pass canceller 
(medium Q) 


2.28 


1 


1 


0.9 


Low pass canceller 
(high Q) 


On each curve, the 


theoretical 


transfer function is 


plotted as 


a continuous 


line . 


The 


data points marked with 


"+" ' s are 


shown also. 


On many of 


the curves , agreement 



between theory and practice is extremely good for the first 
two teeth of the comb. Agreement falls off at higher 
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frequencies for a number of reasons and this will be taken 
up in the next section. However, in a few curves (no. 20, 

21, and 23) the agreement does not seem to be so good. In 
curves 20 and 21, it appears that the coefficient b^^ was 
not set correctly. The correction factor mentioned in the 
first part of this section was used to set all the poten- 
tiometers for these measurements and it is likely that the 
potentiometer itself was not linear throughout its range of 
20 KS^. The data taken seems to be symmetrical about half 
the sampling frequencies and periodic with respect to f . 

So the error in b^^ appears to be the dominant error. In 
curve no. 23, we have the reverse situation ... the actual 
data appears to be taken for a coefficient considerably 
larger than 0.6. Several other curves were taken for these 
levels of feedback with, in general, good agreement. Figure 
2.27 shows an interesting effect for which we have struggled 
to obtain some explanation at least qualitatively. The 
transfer function clearly indicates a rising amplitude 
response with frequency. This rising phenomenon was noticed 
in several curves, but not so pronounced as in no. 27. We 
shall consider several possible sources of this behavior in 
the next section. 

2.2.3 Corrections in the data because of physical instrumentation 
Some effects cause the data to roll off with increasing 
frequency. Other natural effects are periodic and cause dis- 
torted amplitude response. We shall discuss several of these 
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effects which are all caused by the nonideal nature of our 
experiments . 

The most significant cause of rolloff in amplitude 
response of all the first order filters is the sample-and- 
hold effect. Nyquist's sampling theorem states that if a 
signal is sampled at twice its highest frequency component, 
it can be reproduced precisely by passage through an ideal 
low pass filter. We used a sample-and-hold circuit as our 
low pass filter. It is not ideal. The ideal filter ampli- 
tude response is shown in Figure 2.29. On the other hand, 
the amplitude response of the sample-and-hold circuit is 
shown in Figure 2.30. The general shape is sin x/ x where 

X = IT f/f_ (2-25) 

This frequency response must be multiplied by the theoretical 
filter transfer function to get the projected overall trans- 
fer function corrected for sample-and-hold circuitry. This 
has been done for two of the curves shown in Section 2.2.2. 
Figures 2.31 and 2.32 show the result. As can be seen, the 
sample-and-hold factor reduces the theoretical calculations 
substantially (because most of our work is relatively close 
to the Nyquist frequency) and is the major reason for the 
frequency rolloff with increasing frequency. 

The sample-and-hold factor can also be seen intuitively 
by considering the action in the time domain. Our clock 
(sampling) frequency is 20 kHz. We have taken data at 5 kHz, 
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Figure 2.29. Ideal low pass filter. 




Figure 2.30. Transfer function for a sample-and-hold 
circuit, f^ is the sampling frequency. 
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Figure 2.31 Theoretical and Experimental 
a first order recursive filter, coirects 
hold effect, ^ » ^l " * ^1 ” 



transfer function for 
for the sample-and- 
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Figure 2.32. Theoretical and Experimental transfer function for 
a first order recursive filter, corrected for the sample-and- 
hold effect, = 1 , = 1 , = -0.9 
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well below the Nyquist frequency of 10 kHz for which we 
ought to be able to reconstruct a signal. However, even at 
5 kHz, we are taking but four samples per cycle. This 
effect is best seen on an oscilloscope where the four dis- 
tinct levels per cycle can be observed. It is clear that 
the resulting waveform is not a clean 5 kHz sine wave, but 
instead a waveform which contains many different harmonics. 

The second effect which we shall consider is the CCD 

efficiency. As mentioned in Chapter one, as each charge 

packet travels down the CCD, it is degraded slightly because 

of the incomplete transfer at each gate. The measure of 

incompleteness is called inefficiency, e, and is commonly 

-4 

on the order of 10 . The factor for inefficiency is 

exponential and is given by 

(-ne (1-cos 2tt f/f ) ) 

Total inefficiency = e (2-26) 

Using this expression as the attenuation of the signal 
after each passage through the CCD, we can calculate a new 
transfer function for the filter. The block diagram for 
explicit inclusion of this factor is shown in Figure 2.33. 
The overall effect of this factor is shown in Figure 2.34 
for a filter with a^ = 1, a^ = -1, and b^^ = 0.7 (high pass 
integrator) . There are six curves plotted on this figure. 
The curve which is not attenuated at all with increasing 
frequency corresponds to an inefficiency of zero . . . that 
is, the whole charge packet is transferred each time. The 
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Figure 1.33. Explicit inclusion of the factor, F, for inefficiency. 
F = exp(-ne (l-cosZTff/fg )) . First order filter. 
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Figure 2.34, Theoretical transfer function for a first order 
recursive filter corrected for inefficiency, k = ns= , a = 1, 
= -1. = 0.7 ° 
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second curve down from the zero curve corresponds to a value 
of .0125. Each subsequent curve corresponds to an ineffi- 
ciency of .0125 greater than the previous curve. 

There are two important phenomena to notice in Figure 
2.34. One is, of course, the falling off of the maximum 
as frequency is increased. The other, is just as signifi- 
cant, although not so easily seen, and that is the deter- 
ioration of the null (zero) as frequency is increased. This 
is intuitively appealing since, at low frequencies we know 
that a = a, for our filters derived by the bilinear z 
transform. When frequency increases to such an extent 
that the effect of the inefficiency has a considerable 
effect on the gain of the CCD, a^ no longer equals a^^. The 
coefficient a^^ has been shifted by the shift in gain caused 
by the inefficiency. If a^ does not equal a^, then proper 
cancellation at the nulls does not take place. In other 
words, the zero which, at low frequencies, lay at z = -1, 
shifts at higher frequencies so that it lies at z = -z^ 
where |z^l <1. Since the zero is no longer on the unit 
circle, we cannot expect to achieve a perfect zero in our 
transfer function. 

Since the factor for inefficiency is periodic with 
respect to the sampling frequency, we expect that the shape 
of the curve in Figure 2.34 will be reproduced at higher 
regions of the spectrum. Figure 2.35 confirms this expec- 
tation. Calculations for Figure 2.35 have been made from 
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Figure 2.35 Theoretical transfer function for a first 
order recursive filter corrected for inefficiency, extended 
to fg = 20 Khz. k = n , a = 1 , a-j = -1 , = b.7 
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dc to f (20 kHz) . We see the familiar symmetry about 

o 

fg/2, as well as the aforementioned deterioration of the 
zeros and attenuation of the maxima. 

The obvious next step is to combine the effects of 
inefficiency and sample-and-hold circuitry and determine 
what the new theoretical transfer function ought to be. 

This has been done for two representative filters. Figures 
2.36 and 2.37 show different aspects of the transfer function 
for the filter with a^ = 1, a^^ = -1, and b^^ = 0.7. The 
first figure here shows the transfer function from dc to 
10 kHz corrected for both effects. The second figure shows 
the comparison of the experimental data with the new theoreti 
cal transfer function based on an estimated inefficiency of 

-3 

10 . Agreement is good. Figure 2.38 shows the new theoreti 

cal transfer function for the filter with the coefficients 
a^ = 1, a^^ = “1, and bj^ = -0.7. Four teeth are shown. The 
same phenomena are evident with the additional expected 
nonsymmetry about f /16, 3f /16, etc. This nonsymmetry is, 

o o 

of course, caused by the sample-and-hold and inefficiency 
factors, which are not periodic with respect to f_/8 and 
are not symmetrical about f /16. The rising null can also 
be seen with increasing inefficiency. 

One possible explanation for the rising phenomenon men- 
tioned in the preceding section is that the potentiometers 
have an impedance function which falls with frequency. 

These potentiometers are the wire-wound type and one effect 
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Figure 2.36 Theoretical transfer function for first order recursive 
filter, corrected for both sample-and-hold effects and inefficiency. 
Different curves are for different values of inefficiency, as 
explained in the text, a = 1 , a, =-l, b = 0.7 
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Figure 2,37 Theoretical and experimental transfer function for a first order 
recursive filter corrected for both sanple-and-hold effects and an inefficiency of 




Figure 2.38 Theoretical tranrfer function for first order recursive 
filter, corrected for both sanple-and-’nold effects and inefficiency. 
Different curves are for different values of inefficiency as explained 
in the text, a = 1 , a-| = -1 , b = -0.7 
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of such resistors is that they become capacitive at high 
frequencies due to the capacitance between adjacent loops 
of wire. Such a parasitic effect, along with the parasitic 
capacitance of the coaxial cable utilized in this experiment, 
would have the result of reducing the absolute value of the 
feedback coefficient at high frequencies. This, in turn, 
for settings of positive b^^ has the effect of increasing 
the amplitude response. For instance suppose we have 



H(z) 



1 + z 



-1 



1 + 0.3 z 



-1 



The absolute value of H(z) at dc is 2/1.3 or 1.54. At 
higher frequencies, when 0.3 is reduced by a factor of 0.7, 
we have the following transfer function value: 

1 + z"l 

H(z) = ■ Y = 

1 + 0.21 z~^ 

which represents an increase of about 7% which is about what 
is observed in Figure 2.26. For higher values of feedback, 
say 0.7, we have the following two values: 



H(0) = 1.18 H(f^) = 1.34 

where the increase is about 14% which corresponds well 
with the actual observed data in Figure 2.27. 
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strangely, the phenomenon did not occur for = 0.9. 

It is possible that the coaxial connections were changed for 
the data obtained at this high feedback setting. 

To summarize this section, we have discussed the following 
factors which cause deviations of experiments from theory; 

A. Sample-and-hold causes attenuation of 

maxima and nonsymmetry of comb filter 

B. Inefficiency causes attenuation of maxima 

and deterioration of zero (null) . Also causes nonsymmetry 
of comb, but its effect on symmetry is not so marked as 
the sample-and-hold effect. 

C. Potentiometer causes an increasing of 

amplitude response for positive feedback cases with increasing 
frequency. Causes a decreasing of amplitude response for 
negative feedback case with increasing frequency. 

2.2.4 Selectivity 

The idea of selectivity is central to filtering. A 
filter changes the frequency content of a signal. Highly 
selective filters eliminate almost all frequency components 
except the small band which they pass. No discussion of 
filters is complete without some mention of selectivity. 

The basic measure of selectivity is Q. Q is defined in 
various places by various equations, but it is appropriate 
here to use the definition 



^center 

bandwidth 



(2-27) 
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Figure 2.39. Frequency response of potentiometers. 



where refers to the center frequency of the passband 

and the bandwidth is the distance (in frequency) between 
the two 3 db down points. 

The above definition is fine for integrators. The 
center frequency is well defined and the bandwidth is small 
(usually) relative to the center frequency. The above 
definition breaks down when applied to cancellers. In 
canceller type filters, the interest is not really in the 
passband, but in the stopband. A known noise signal exists, 
is received, and filtered. The canceller ideally stops all 
the noise. Of course, in practice, this is not possible 
for some of the noise components are in the passband of the 
filter. Nevertheless, if we have a good idea of the spec- 
trum of the noise, then we can build a filter which is opti- 
mum for the noise under consideration. At first, we attempted 
to use the reciprocal of the definition given in Equation 
(2-27) . This did not work so well since it made Q a non- 
linear function of the center frequency. A better choice 
of Q to be used with cancellers appears below: 

A. H(f) is the transfer function of a canceller 

B. Normalize H(f) to unity 

C. Form the difference 1 - H(f) = H'(f) 

D. Now utilize Equation (2-27) and find the Q of 
H* (f) . 

This is the Q of H(f) as well. 

This procedure retains the linearity and has worked well 
in our calculations. 
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Figure 2.40 shows the theoretical calculation of Q for 
a high pass filter derived by bilinear z transform with 
variable feedback coefficient. As can be seen, the Q is 
low for negative values of b^^ and also for positive values 
of bj^ less than about 0.7. Q rises dramatically for high 
values of b^^ (near 1) . Of course, the usual definition of 
Q was used in this calculation. The use of the 1 - H(f) 
definition would produce the identical curve except that 
the rise in Q would occ\ar near -1. The overall evaluation 
of a design should use the appropriate definition of Q. 

2.2.5 Summary 

In this section, we have presented our experimental 
results together with theoretical calculations and some 
reasons why the experimental data do not always agree 
exactly with the theory. The first order filter has the 
advantage of simplicity and the use of the bilinear z trans- 
form often provides us with a rather large maximum- to-minimiim 
response ratio. One of our conclusions is especially 
significant; The standard z transform is a rather poor 
tool to use when designing filters or other processors with 
a frequency response specification in mind. Since most 
design is done to frequency specifications, we would strongly 
suggest that the standard z transform not be chosen for a 
means to design unless the sampling frequency is much much 
greater than the sampled analog cutoff frequency desired. 
Several factors interfere with perfect operation of our filters 
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Figure 2.40 Quality factor (Q) as a function of feedback 
coefficient for a high pass first order recursive filter. 
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A. Instrumentation problems such as the frequency 
response variation of the potentiometers. 

B. Inefficiency of the CCD. 

C. Nonideal response of the sample-and-hold 
circuitry. 

Figure 2.41 and 2.42 are of interest. They summarize 
the theory of first order sampled analog filters designed 
with the bilinear z transform. Both figures show the same 
thing, but from different aspects. The figures show the wide 
range of transfer functions available to a designer merely 
by changing the feedback coefficient, b^^, of a first order 
filter. The x-axis is the b^ value, the y-axis is frequency, 
and the z-axis is, of course, the amplitude response. The 
coefficient b^ varies from -0.9 to 0.9 while frequency varies 
from dc to 2.5f^. Where the curves have especially wide 
passband is where b^ = +0.9 (foreground in Figure 2.41, 
background in Figure 2.42). The other extreme (b^^ = -0.9) 
is characterized by a very narrow passband and wide stopband. 

SECTION 2.3 THE RECURSIVE SECOND ORDER COMB FILTER 

2 . 3 .*1 The basic characteristics of the second order comb 
filter are similar to those of the first order except that 
two poles and as many as two zeros may occur in the transfer 
function. The sampled analog transfer function is 

u- -1 u- -2 

a + a, z + a.,z 

H(z) = -2 £_ (2-28) 

1 + b^z + b2Z 
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Figure 2.41. Low pass theoretical transfer function for a first order 

recursive filter with various values of feedback, a = 1 , a = 1 

0 1 
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Figure 2.42 Low Pass theoretical transfer functions for a first 
order recursive filter with various values of feedback, a = 1 

a, . 1 
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The second order recursive filter is widely used for such 
applications as radar signal processing, bandpass and band- 
stop filtering, and in coinmunications . It is more widely 
used than the first order filter since it offers a wider 
variety of transfer function shapes for a given sampling 
frequency . 

The comb filter design technique, as in the first order 
case, utilizes a serial-in serial-out delay line, such as 
the CCD discussed earlier. The true transfer function for 
our experiments using two eight stage delay devices is 



H(z) 



a + aoZ ^ + a,-z 
o 8 16 

1 + bjz'® + 



(2-29) 



Such a transfer function is clearly a sixteenth order filter, 

one with, at most, sixteen poles and sixteen zeros. As 

before, there are many coefficients missing; a^,a 2 ,...f 

a^,ag, . . . . ,a ^2 and the same corresponding b's. This design 

procedure enables the engineer to consider his sixteenth 

order filter as a second order filter where his adjusted 

sampling frequency is one eighth that of the true sampling 

frequency. As shown in Appendix C, the poles and zeros 

for such a transfer function are not equidistant as in the 

first order case. Figure 2.43 shows that the poles and zeros 

are arranged so that they are symmetrical with respect to 

dc, f_/8, 2f_/8, nf_/8. Such a symmetric nature 

s s s 

yields a transfer function periodic with respect to fg/8. 
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Figure 2.43. Symmetrical pole plot on the z-plane. 
Poles are spaced symmetrically to radial lines of 
the form nf^/8, where n=0,l ,2,. .. ,7. This affords 
the basis by which we make the sixteen- to-second 
order translation. 
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The ability to regard a sixteenth order filter as a 
second order filter is important and significant for the 
same reasons which applied to the eighth-to-f irst order 
translation. However, a second order filter is not so 
simple as a first order filter. The effect on the transfer 
function of a variation of one or more coefficients is not 
generally intuitively obvious. For instance, an engineer 
may need to know the effect of the variation of b^^ and b 2 
by some small percentage because of thermal variations. 

Such parameters as cutoff frequency, "skirt slope", and 
ripple must remain within system tolerances, even though 
the coefficients vary. An analysis of these effects has 
been made and will be explained fully in Section 2.4. The 
analysis is for second order recursive sampled analog filters. 
The same analysis undertaken for a sixteenth order filter 
would, of course, be a huge project. 

Section 2.4 DESIGN METHODS FOR A SECOND ORDER SAMPLED 
ANALOG FILTER 

2.4.1 Introduction 

The design of second order filters is necessarily more 
complex than the design of first order filters. The first 
order filter had but two possible singularities, one pole 
and one zero, and, moreover, each was by necessity, real. 

Thus , we were dealing with singularities restricted to the 
real line on the z-plane, and design was essentially a 
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one-dimensional problem . . . deciding where on the real line 
to place the singularities. Obvious rules were available ... 
as shown in the experimental and theoretical figures of Sec- 
tion 2.2/ the effects on the shape of the transfer function 
by the position of the pole and zero become intuitive and 
easy to work with for the design engineer. However/ the 
relationship between the shape of the transfer function and 
the position of the poles and zeros of a second order filter 
is not so well understood. Of course, as the order of the 
filter increases, the effect of the placement of the poles 
and zeros becomes more and more complex. 

Therefore, the design procedure for second order and 
above can no longer rely on intuitive means. Analytical 
processes are required for the proper design of second order 
filters. There are to be discussed here three different 
methods of designing such filters: 

A. Replication in the z-domain of an analog filter 
whose characteristics in the s-domain are well known and 
well documented. 

B. Selection of coefficients by reference to charts 
which indicate the positions of the relevant parameters: 
cutoff frequency, ripple values, etc. In this section we 
will derive and discuss the b-plane triangle of stability 
for second order sampled analog filters. 

C. Direct sampled analog filter design using cost 
functions and the general purpose digital computer (Fletcher- 
Powell algorithm) . 
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Section 2.4.2 Replication in the z-domain of known analog filters 



Three of the best known analog filters are the following: 
Butterworth, Chebyshev, and elliptical (sometimes called 
Cauer parameter filters) . Briefly, the Butterworth filter 
has a frequency characteristic which is monotonic in both 
passband and stopband, has steeper skirts for higher orders, 
and has a maximum error (deviation from ideal) in each stop- 
band and passband which occurs at one and only one frequency. 
This maximum error gave rise to the thought that it might be 
possible to "spread" this error throughout a range of frequen' 
cies, and thus make the error at any one frequency much less. 
This was accomplished by the use of polynomials in the denom- 
inator of the transfer function which were called Chebyshev 
polynomials. Polynomials which "spread" the error throughout 
the passband only have roots (poles of the transfer function) 
which lie on an ellipse in the s-plane. This location of 
poles results in a designed-in ripple effect in the passband. 
The size of the ripple is one of the parameters of the filter 
and is held constant by a specification, usually called e. 
This specification often determines the order of the filter 
required since higher orders can reduce the size of the 
ripple. 

The elliptic filter takes advantage of these polynomials 
to produce a ripple in the stopband as well. In other words, 
it spreads the overall error in the stopband over a wide 
range of frequencies. 
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2. 4. 2.1 The Butterworth Filter 



The Butterworth filter is characterized by having a 
monotonic transfer function in both the passband and stop- 
band. The squared magnitude transfer function is given by 

|F(w)|^ = (2-30) 

1 + (JiL)2n 
‘^c 

where n is the order of the filter. The design procedure 
is straight-forward. The poles are located on a circle of 
radius w and are equally spaced. If the order is odd, 
then one pole lies on the negative real axis. If the order 
is even, then the negative real axis lies midway (angularly) 
between two poles. 

For first order, there is one pole and no zeros. The 
transfer function was considered in Section 2.1. Both the 
standard z transform and the bilinear z transform were applied 
to it to derive a replication in the z-domain. As was noted 
at the time, the standard z transform produced impulse invari- 
ance while the bilinear z transform insured that the ranges 
of the two transfer functions (the analog and the sampled- 
analog) would be precisely equal over the domain of real 
frequencies. Warning was given for the standard z transform 
because a designer could not be sure at the outset whether the 
cutoff frequency, known for the analog filter, would even 
exist in the sampled analog filter. Expressions for the 
cutoff frequency were given in terms of the coefficient, 
b^f and inversely. 
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The second order Butterworth filter can be transformed 
in exactly the same ways . . . either by standard z transform 
or the bilinear z transform. The second order Butterworth 
polynomial has a transfer function of the form 

H(s) = j (2-31) 

s +/2 sto + 0 ) 
c c 

The poles of the transfer function lie at the following 
positions: 



s 





(2-32) 



Once these things are known, the design of the corresponding 
sampled analog filter can proceed. First we consider design 
by the impulse invariant method. The first step is to find 
the residues and form 



H(s) = Z 



A. 

X 



. , s + s . 

X=1 1 



(2-33) 



Then the inverse LaPlace transform is found as before 



h(t) = Z A. e 
i=l ^ 



-s . t 
1 



(2-34) 



This expression is sampled at some interval, T, and the 
Z transform taken of the resulting sampled signal: 
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Z(h* (n) ) = I 



(2-35) 



1 



e sin 3T z~^ 



— otT — 1 

e 2cos 3T z + 



-2aT 

e 



-2 

z 



which means that the coefficients of the canonical form for 
a sampled analog filter are 



“ 

a^^ = ^ e sin 3T 
^2 = 0 



b 

o 




= 1 

— OtT 

= -2e cos 3T 

_ -2aT 
= e 



In order to find the cutoff frequency of the analog transfer 
function, we form the equation 



|h(w^) = |h(0) 1^/2 



(2-36) 



These equations must be solved. The results are given in 
Equation (2-37) ; 



f 



2 

c 





(2-37) 



where the factor co^ is the same quantity as in the Butterworth 
polynomial. The roots of this expression are four in niomber 
including one positive real root, one negative real, one 
positive imaginary, and one negative imaginary. We naturally 
take the positive real root as our cutoff frequency. 

For the sampled analog case, the cutoff frequency Ccinnot 
be calculated so easily. The form for the equation relating 
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the cutoff frequency to the coefficients can be derived and 
solved but it is tedious. The expression is a quadratic 
and includes the eigenvalues of the analog transfer function 
as parameters: 



: _ 1 cos gT + cos gTe'^^*^ ± Jcos ^ BT ( l+2e"‘^“'^+e"®“'^) 

c 2 ttT ^ 



• + cos BT[-8e““'^ + -8e"^“'^]+l+10e"^°‘'^+e“‘^“'^ } 

(2-38) 

This expression, though complicated, is straightforward. 
Simply insert the values of a and g which are the coordinates 
of the poles and zeros in the s-plane and compute. That 
these values are always known is inherent in the ability to 
design an analog Butterworth filter. Naturally, the solution 
to this equation can perhaps be done best by a computer. A 
program for its solution was written and tested. 

In order to use the bilinear transformation, we need not 
find the residues or find the inverse LaPlace transform of 
the analog transfer function. We make the substitution as 
we did in Section 2.1 for a first order filter: 

2 z - 1 
® T z + 1* 

The second order Butterworth filter, transformed , via the 
bilinear z transform into the z-domain produces the following 
results: 
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H(z) 



(2-39) 



1 + 



G(z~^ + 2z ^ + 1) 

-4+2a^+26^, _-l ^ ,4-4a+a^+B^^ _-2 

— 2 — 1' ^ ^ y ' ■> ' ^ 

4+4a+a +3 4+4a+a +3 



where G is the overall gain of the filter and is given by 

G = 2- (2-40) 

4 + 4a + a +3 

The coefficients of the sampled analog filter are easily 
picked off the transfer function: 



a 



o 



1 



a 



1 



2 



b = 
o 



1 



-4+2a^+23^ 
2 2 

4+4a+a +3 



^2 



1 



4-4a+a +3 
2 2 

4+4a+a +3 



where a and 3 are the real part and the imaginary part of 
the pole pair (s-plane) respectively. 

The cutoff frequency can be fovind in the usual manner. 
Like the result in the impulse invariant approach, the cutoff 
frequency calculation turns out to be complicated. It is a 
function of the singularities in the s-domain. If we let 
the cutoff frequency be a f\anction of the sampled analog 
coefficients then we obtain a quadratic equation in 
cos (2tt f /f ) for which the following are coefficients: 

C o 
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term coefficient 



2 

cos 


(2it 




1 


+ + 


2bib2 + 


2b^ - 


6b2 


cos 


(27T 


V^s> 


2 




+ 4bj 




(2-41) 


o 

cos 


(2tt 


V«s> 


-1 


2 2 
- b^2 . ^ 


2b,b2 + 


2bj^ + 


2*=2 



Notice that these equations are perfectly general in 
that they don't appear to require the use of the Butterworth 
model to be correct. As a matter of fact, we will find that 
these same equations apply in the cases of the Chebyshev 
and elliptic filters, with the difference being that the 
coefficients are derived in different ways. 

2. 4. 2. 2 The Chebyshev filter 

The design of a Chebyshev filter requires the specifica- 
tion of the maximum ripple in the passband, and the minimum 
attenuation at some frequency greater than the cutoff fre- 
quency. The order of the filter will be derived from these 
considerations. Once the order is found, the poles of the 
transfer function are found from the major and minor axes 
of the ellipse: 

The -radius of the minor axis is aw^ where 
1 _1 

a = i (y^ - Y Y ~ ^ (2-42) 

and the radius of the major axis is b^^^ where 
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(2-43) 



1 _1 
= j (Y^ + Y 

Thus a and b are functions of known quantities, the maximum 
permissible ripple given by e, and the order of the filter, 

n. 

Since there is no ripple in first order filters, it makes 
little sense to talk about first order Chebyshev filters. 
Nevertheless, the transfer function for such a filter would 
be 



H(s) = i (2-44) 

1 + es/s^ 

where the only pole is clearly on the real axis of the s-plane. 
There are no zeros, which is the same condition as exists 
in the Butterworth filter family. 

The cutoff frequency for the transfer function in Equation 
(2-44) is 



f 

c 



!i 

2tt £ 



(2-45) 



Now if we were to calculate the transfer function for the 
sampled analog filter based on this Chebyshev model by the 
standard z transform approach, we would make the following 
calculations : 

A. Find the residues, as was done for the Butterworth 

filter. 
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B. Take the inverse LaPlace transform 

C. Sample the impulse response 

D. Take the Z transform of the sampled impulse 

response. 

The result for first order is 



H(z) 




(2-46) 



9 

Given in this form, it looks like what was encountered in 
Section 2.1, except in this case, the exponential in the 
denominator has, in its argument, the maxim\am ripple factor 

e . 

The sampled analog cutoff frequency turns out to be 



f 

c 



1 

2ttT 



cos 



-1 




4bj_ + 



2b, 



(2-47) 



which is the same equation as we had in Section 2.1. 

The same Chebyshev filter could be transformed by the 
bilinear z method with a resulting different set of sampled 
analog coefficients. Making the usual sxabstitution 



2 z - 
T z + 1 



we get a transfer function which is very similar to the 
transfer function of Equation 2.4 except that the maximum 
ripple factor, e, is once again included in the denominator. 
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The coefficients are as follows; 



a 



o 



1 




1 



^1 = 



1 



s . T - 2e 
1 

s .T + 2e 
1 



Because of the requirement for stability, the coefficient, 
b^, is limited to absolute values between zero and one. For 
these values, the cutoff frequency can be found. The equation 
for the sampled analog cutoff frequency is given in Equation 
(2-5) in Section 2.1. It is repeated here for clarity: 



^c = 



1 

2ttT 



cos 



-1 



-2b, 



1 + b. 



(2-5) 



As before, this equation can be solved for b^^ in terms of 
^c 



-1 ± sin 2irf T 

b = (2-6) 

cos 2irf^T 



The second order Chebyshev filter does exhibit ripple. 

The passband contains a ripple no greater than e, and the 
stopband is monotonic. The second order Chebyshev polynomial 
is 



P2(s) = 2s^ - 1 



(2-43) 
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The form of the squared magnitude of the filter is 



|h,(s) 



2 2 

1 + e P (s) 
n 



(2-49) 



where = 1, = s, ?2 = 2s - 1 and, in general, 

= 2sP^ - Often P is given in terms of s/s^. 

In order to design a second order Chebyshev sampled 
analog filter, it is necessary, as always, to find the 
inverse LaPlace transform of the H 2 (s) transfer function 
(for the impulse invariant method) . This means finding the 
residues at the two poles of the transfer function. The 
poles are found in the same way as in the first order with 
the additional complexity being that the poles are complex, 
not real. Each pole lies on the ellipse defined earlier 
in this section. The position of the poles on the ellipse 
is found by dividing the circle of revolution in the s-plane 
into 2n angular sectors, each with angular measure fr/n. If 
n is odd, then one boundary of the sectors lies on the positive 
real axis. If n is even, then the positive real axis lies 
directly in the center of one of the sectors. This is illus- 
trated in Figure 2.44. 

The precise position of the poles are then found by the 
following procedure; 

A. The abscissa of the coordinates of the pole is 
specified by the point on the circle defined by the minor 
axis of the ellipse, and through which the boundary of the 
sector passes. 
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Figure 2.44. Location of poles for a third 
order Chebyshev filter. 
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B. The ordinate of the coordinates of the pole is 
specified by the point on the circle defined by the major 
axis of the ellipse, and through which the boundary of the 
sector passes . 

After the poles are found, the transfer function must 
be arranged in the residue form, inverse LaPlace transform 
taken, sampled, and the Z transform taken. This procedure 
(now standard) results in the following transfer function: 



Z(h(n)) = 



2/B 



-aT 

e sin 



6T z 



-1 



- 2e '^'^cos BT z ^ 



~ -2aT -2 
+ e z 



(2-49) 



and the sampled analog coefficients can be easily seen from 
the above form: 



a = 0 = 1 

o o 

a^ = 2/B e sin BT b^ = -2e cos BT 




where a and B are the real and imaginary parts of the 
singularities in the s-plane. 

The cutoff frequency for this analog filter can be found 
by determining for which frequency the squared magnitude 
function is down by 3 db. This entails solving the following 
equation: 
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1)^ = 2(1 + e^) 



(2-50) 



2 2 
1 + e^(2(;^)^ - 

“i 

which, of course, yields four roots. For any real e, the 
equation will yield one positive root, one negative root, 
and two purely imaginary roots. Naturally, it is the positive 
real root which we want for our cutoff frequency. For 
e = 0.01, for instance, the squared roots of the equation are 
found by 



& 2 = (2-51) 

00 . « 

X 2 

which yields roots at +50.5 and -49.5. The square roots of 
these numbers are approximately ± 7.1 and ±j7.0. Thus the 
cutoff frequency is about 7.1 

For the sampled analog case, the cutoff frequency can 
be found in the same manner. The resulting expression is 
identical to that found in the section on Butterworth filters: 

, _ 1 -l,cos BT ( l+e'2“’') i ) cos 26T (l+2e“‘“’'+e-^“’’) . . . 

^c 27 tT ^ - -2aT 

2e 



+ cos6T[-8e~°'^ -f 4e~2°‘'^ - 8e~2°‘'^].|.H-10e~2°‘'^te~^°‘'^ ) 



2e 



-2aT 



(2-52) 



The quantities a and g are, of course, different for the 
Butterworth and Chebyshev cases. 
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In order to use the bilinear z transformation, we need 
not determine residues or find the inverse LaPlace transform 
of the transfer function. We make the usual substitution 
and find that the second order Chebyshev low pass filter in 
the z-plane has the following transfer function: 



H(z) 



G(l+2z~^+z~^) 

T ^ (-4+2a^+26^) -1 . (4-4a+a^+6^) -2 

1 + = = — z + 2 — 2 — z 

4+4a+a +6 4+4a+a +6 



and the coefficients are 



(2-53) 



% = 1 



b = 1 
o 



ai = 2 



b, = 



_ -4+2g^+2B^ 



4+4a+a^+3^ 



(2-54) 



^2 ^ 



b^ = 



4-4a+g^+e^ 
4+4g+g^+3 ^ 



as found before in the Butterworth section. 

The cutoff frequency is found in the usual manner. The 
same expression is obtained as in the Butterworth section. 
The expression is a quadratic in cos {2ir f_/f„) with the 
following coefficients; 



term 



cos {2ir f /f ) 

w o 



cos {2ir f /f ) 
c s 

cos°(27t 



coefficient 

1 + + 2bj^ + 2b^b2 - 6b2 



2 + 2h^ + 2b2^ + 4b2 



(2-55) 



-1 - - b2^ + 2b^b2 + 2b^^ + 6b2 
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Once again this is a complicated expression, but given 
the poles, it is a straightforward procedure to find, first 



the coefficients, and then f . 

c 

2. 4. 2, 3 The Elliptic Filter 

The elliptic filter (Cauer parameter filter) is a natural 
outgrowth of the Chebyshev filter. The Chebyshev filter is 
characterized by having ripples in the passband and a mono- 
tonic characteristic in the stopband. The elliptic filter 
"spreads" the error in the stopband of a Chebyshev filter by 
producing ripples in the stopband as well as in the passband. 
The elliptic filter utilizing this procedure produces the 
best (shortest) transition region for any filter of a given 
order and given ripple amplitude. In other words, the 
elliptic filter produces the steepest skirts from passband 
to stopband. 

The specifications of an elliptic filter must begin with 
the maximum ripple in the passband which is defined by e. 

The maximum ripple in the stopband is given by 1/A. In 
addition, the cutoff frequency, f , and the cutin frequency, 
f^, must be specified. 

The squared magnitude form of the filter is given by 




1 



(2-55) 



2 2 

1 + £ F^(S) 



where F is an elliptic function. 
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The elliptic function, sometimes written as sn(u,k) 
must be solved to find each pole and zero location. The 
expression sn(u,k) is shown below where u is the independent 
variable and k is a parameter. If y = sn(u,k) then 

y 

u = / — = sn"^y (2-56) 

° (1-t^) ^/^(1-k^t^) 

The form of F can be understood by the following two forms : 

If the filter is to be of even order, then 

, 2 2, , 2 2, , 2 2, 

(tO^ -W ) (C0_ -W ) ... (cO- ^ -CO ) 

F(to) = - 7 o T- (2-57) 

(l-co^co^^) (l-co^t02^) ... (1-w “2n-l ^ 

If the filter is to be of odd order, then 

, 2 2 , , 2 2 , , 2 2 , 

(co^ -CO ) (co. -CO ) ... (co, -CO ) 

F(co) = 2 — ^ (2-58) 

(I-CO 2 CO ) (1-co^^co ) ... (l"‘*^2n ^ ^ 

The trick is to choose the zeros of F (and thus, the 
poles since if there is a zero at co^, then there is a pole 
at 1/co^) such that the maximum ripple e, the minimum atten- 
uation A, and the f and f are fulfilled. To do this requires 
solving the incomplete elliptic integral shown in Equation 
(2-56) . 

In order to calculate the coefficients of a sampled 
analog filter, by the impulse invariant method, it is first 
necessary to find the poles of the analog transfer function. 

The poles must be found by solving 
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1 + e^sn^(u,k) 



0 



(2-59) 



No rigorous attempt has been made to solve this equation. 
Values for various arguments of sn(u,k) have been calculated 
and are available in tables. 

After the poles are found, the derivation of the seimpled 
analog filter types continues along the same lines as for 
the Butterworth and Chebyshev filters with the single excep- 
tion that the analog filter form contains zeros. This fact 
does not appear to be significant, however, for purposes of 
design in the z-plane. 

2.4.3 DIRECT SELECTION OF COEFFICIENTS FOR THE SAMPLED 
ANALOG SECOND ORDER RECURSIVE FILTER 

The derivation of a Scimpled-analog filter from well known 
analog filter types is attractive because some of the quali- 
ties of the analog filters are attractive. Unfortunately, 
upon transformation, many of their attractive traits are 
lost. Certainly, by use of the standard z transform, one 
does not even have the assurance that the ranges of the 
analog filter and the sampled analog filter are the same. 

On the other hand, the frequency "warping” effect spoken of 
in Section 2.1 which occurs when the bilinear z transformation 
is used often is a hindrance to proper design. 

It appeared, in our research, that what was really needed 
by the design engineer was some method to determine, first, 
the admissible values of, say, b^^ and b 2 in a sampled analog 
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filter; second , the general shape of the transfer function 



as a function of and h2t and third, some indication as 
to how much critical parameters change when the coefficients 
are shifted slightly due to temperature variations, aging, 
etc. This last consideration is especially important since 
the recursive filter has been held in some disfavor by 
designers because of possible stability problems, alleged 
problems of xanstable parameters even if the filter itself is 
stable, and finally, because of its infinite duration unit 
pulse response (time domain) . If it could be shown that 
the filter of a certain design was stable and had stable 
parameters over a specified tolerance of coefficients, then 
the idea of a recursive filter in systems should be encouraged. 

To that end, an effort was made to discover the admissible 
values of b^^ and b2 in a second order sampled analog recursive 
filter, whose transfer function is as follows; 



There are three criteria defining admissibility of a point 
whose coordinates are (b2/bj^). One criterion is the stability 
constraint, a physical constraint. Thus, if is a pole of 
the system, then 



H(z) = • 




-2 



-2 



( 2 - 59 ) 




( 2 - 60 ) 
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The other two criteria defining admissibility arise from 
basic algebra. For the quadratic equation 




where and are roots to the equation / then 




(2-61) 



and 




b 



1 



(2-62) 



For purposes of clarity, it is appropriate to treat the 
admissibility problem in three distinct cases. For the 
first case, we set b 2 negative. 



If b 2 is negative, then the polarity of the poles (their 

signs) are opposite. Using the third criterion (Equation 

(2-62)), we observe that the absolute value of b^^ must be 

less than the larger of the absolute values of and Z 2 * 

A more useful observation for this case is that there can 

never be complex roots to the polynomial. The discriminant 

2 

of a quadratic equation, ax + bx + c = 0 is 



A. CASE ONE: b 2 negative 



b 



2 



4ac 
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I 



and, in our case, the discriminant is 




Thus, if b 2 is negative, then the discriminant is positive 
for all values of both b^ and b 2 , thus insuring that each 
root is real. 

A graph can clearly show the first two criteria: Figure 

2.45 shows the and Z 2 axes with a hyperbola 

z^Z 2 = b 2 (b 2 is negative) (2-63) 

Also superimposed onto the graph is the xinit square. Since 
neither pole can have an absolute value greater than 1, the 
poles must lie within the square. Combining those two cri- 
teria, we have the darkened portion of the hyperbola drawn 
in Figure 2.45 which is the locus of the only available pole 
locations. Figure 2.46 shows the admissible values of b^ 
based on the data in Figure 2.45. (If z^^ = 1, then Z 2 = ~b 2 
and b^ = -(l-b 2 ) etc.) 

Figure 2.47 shows the locus of admissible pole locations 
on the z-plane for this case. Note that one pole varies from 
b 2 to 1 while the other varies from -1 to “b 2 . Finally, 

Figure 2.48 shows the important result; the permissible 
values of b^^ with respect to b 2 for b 2 negative. The diagonal 
lines have the following equations: 
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(upper line) 



(2-64) 



b 



1 



1 + b 



2 



b 



1 



-1 - b 



2 



(lower line) 



(2-65) 



B. CASE TWO; 



b 2 is positive and the poles are 



real. 



We initially note that the relationship between the poles 
and bj^ is 



This time, since b 2 is positive, each pole must have the same 
sign. Thus the absolute value of b^^ is always greater than 
either z-^ or Z 2 * We also note that, in this case, the dis- 
criminant is not necessarily positive. Nevertheless we 
shall enforce this condition artificially and take up complex 
roots in CASE THREE. 

The hyperbola is once again shown in Figure 2.49. The 
darkened segments of the hyperbola indicate the admissible 
values for the poles given some positive b 2 • Figure 2.50 
shows the admissible values of b^^ based on the admissible 
locations of the poles. 

Going further. Figure 2.51 shows the admissible pole 
locations on the z-plane and Figure 2.52 extends Figure 2.48 
to show the admissible values of b^ given a positive b 2 
and real roots. 




b 



1 



(2-62) 
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The region is bounded by the extension of the two lines 
given in Equations (2-64) and (2-65) plus the horizontal 
parabola 



b^ = 2 /EJ (2-66) 

C. CASE THREE; b 2 is positive and the poles are 

complex. 

For complex poles, the job is made easier by the obser- 
vation that b 2 is always the square of the absolute value of 
each of the two complex conjugate poles; If z^^ = x + jy 
and Z 2 = X - jy, then the absolute value of either of the 
poles is 



/ 2 ^ 2 
7x + y 

Applying the algebraic criteria, we have 

b^ = - 2x 

2 2 
b2 = x*^ + y*^ 

which, as stated, is the square of the absolute value of 
either pole. Since we know that the absolute value of each 
pole cannot exceed one, we can limit our search to values of 
b 2 less than 1. 
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Figure 2.53 shows a root locus diagram (same as Figure 
2.51 except extended to include complex roots). As the pole 
locations are moved around in the b 2 circle, changes con- 
siderably. As a final graph, Figure 2.54 shows the total 
admissible region for b 2 and b^ including all three cases 
discussed here. Notice that the region is convex (the complex 
case having filled in the formerly concave region) and is 
bounded by three line segments: 



The full meaning of Figure 2.54 can be realized if an engineer 
is visualized in the laboratory attempting to design some 
second order recursive filter. He cannot choose his b^^ and 
b 2 at random. He must choose them so that the point (b 2 /b^) 
lies within the region shown in Figure 2.54. We have labeled 
this region the b-PLANE TRIANGLE OF STABILITY since at all 
points on and outside the triangle, the filter becomes 
unstable. 

It remains to be determined how the triangle itself can 
be partitioned. Certianly, the triangle gives us the region 
of real roots and of complex roots, but their significance 
cannot be grasped completely without further partitioning. 



b 



1 



1 + b 



2 



(2-64) 




(2-65) 



b 



2 



1 



(2-67) 
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Figure 2.53. Loci of admissible pole 
locati ons 
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It was with this in mind that we began to investigate a 
filter's characteristic based on a region in the b-plane 
rather than on the pole-zero plot on the z-plane. For this 
test, we held the coefficients a^, a^^, and B 2 at the following 
values: (those naturally arising from the bilinear z trans- 

form) 



^o = 1 




These values of a's insure that zeros will exist in the 
transfer function of the sampled analog filter. This is 
important for the radar processing application we had in 
mind and which is discussed in greater detail in Chapter 
three. The condition a^^ = +2 indicates a low pass filter 
(second order zero at z = -1) and the condition a^^ = -2 
indicates a high pass filter (second order zero at z = 1) . 

A computer program was written to calculate the cutoff 
frequency for many different combinations of b 2 and b^^, 
insuring that the stability criterion was met. This progreim 
is labeled STACY15 and is listed in Appendix D. Figure 
2.55 contains the results of these programs. Figure 2.55 
shows cutoff frequency as a fraction of sampling frequency 
for low pass filters, while Figure 2.56 shows the cutoff 
frequency as a fraction of sampling frequency for high pass 
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Figure 2.55. Isoccntours for cutoff frequency of a second 
order lew pass filter. 3^=1, a^=2, 9^=1. The b-plane 
triangle of stability. 
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-1.5 



Fiouro 2.55, Isccontours for cutoff frequency of a second 
order high pass filter, a =l,a =-2, a„=l. The b-plane 
triangle of stability. o 1 2 
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filters. The triangle of stability has been divided into 
ten sections. Each section consists of b-pairs such that 
the cutoff frequencies lie in a certain range. For instance, 
for extremely narrow low pass filters, we would want our 
coordinates to lie in the crosshatched region of Figure 2.55, 
for this is where the cutoff frequency is less than 0.05fg. 

Each solid line is an "iso-cutoff frequency line" ... a 
line whose locus of points defines filters with equal cutoff 
frequencies. Of particular interest is the solid line whose 
cutoff frequency is 0.25f . This is the dividing line between 
integrators (below the line) and cancellers (above the line) . 

As can be seen from the triangle of stability, large 
latitude for the selection of b^^ and b 2 exists when seeking 
a narrow band low pass filter. On the other hand, only a 
small region exists for a canceller of equivalent Q (cutoff 
frequency greater than 0.45f ). 

Several other observations may be made: 

A. When complex roots are used, there is, in 
general, a maximum of the transfer function somewhere be- 
tween dc and f /2. When real roots are used, the maximum 
s 

exists at dc (for low pass filters) or f„/2 (for high pass 
filters) . The consequence of this result is that if we 
require our filter to be monotonic in the passband, then, 
generally, we require real roots. Such a conclusion is not 
completely valid, since there is a small region where even 
with complex roots, the transfer function is monotonic in 
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both passband and stopband. Now real roots work pretty 
well for integrator type second order filters, but less 
well for canceller types. The narrow region between the 
complex-real roots dividing line and the boundary of the 
triangle for canceller type filters probably precludes the 
design of a useful second order canceller (low pass) using 
real roots. The integrator type has no such limitations 
as can be seen on Figure 2.56. 

B. Rough estimates of how the cutoff frequency 
varies with b^^ and b 2 can also be seen from Figures 2.55 
and 2.56. Consider Figure 2.55 at about = 0.15. The 
cutoff frequency is O.Sf^, and hardly changes at all for 
values of b^^ from 0.5 to 1.2. Thus we would say that the 
directional derivative in the b^^ direction for the cutoff 
frequency in this region is zero. 

All the contours have these regions of nearly constant 
cutoff frequency. 

C. Some contours show marked insensitivity to 
changes in b 2 as well. Consider the 0.25f^ contour. For 
bj^ = -0.4, it is shown that the cutoff frequency hardly 
changes for values of b 2 from 0.5 to 0.75. 

D. The direction of increasing f for Figure 2.55 

is toward the upper right hand corner of the b-plane triangle. 

Figure 2.57 (ripple) shows the dependence of ripple on 
values of (b 2 fbj^) . As can be seen, the entire region of the 
triangle where real roots are found is in a region of zero 
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r = 45% 



Figure 2,57. Isoccntours for ripple (r) of a second 
order low pass filtt^r. a ^ 2^=2, a 2 =l. The b-plane 
triangle of stabiliv.y ^ 
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ripple. Moreover, the zero ripple region extends into the 
complex root region as shown in Figure 2.57 (the region 
between the complex-real root contour and the zero ripple 
contour) . For this calculation, ripple is defined by 



r 






( 2 - 68 ) 



Such a definition permits, at times, certain undesirable 
mathematical results. For instance, for exceptionally 
large ripple, it is possible that the amplitude response 
at the cutoff frequency, lH(f ) |, is actually larger than 
the response at zero, 1 h( 0) ]. To see how this arises, we 
need to have a precise definition of the cutoff frequency 
used in this research. That definition is 




.707 |H (f) 1 + lH(0) 1 
I max ' I ^ ' I 

2 



(2-69) 



which uses, as a basis, an average value of the amplitude 
response in the passband. Now whenever ]H(f ) | = lH(0) | 
we have 



2lH(f^) 1 = .707 (iHj^^^(f) 1 + 1H(0) \) 



or 



1.293 I = .707 | 
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and 



|H(0)| = .547 



SO that the ripple is 



r = 






max 



= .227 



2 l»max<f> 



(2-70) 



Low pass filters with a ripple of 22.7% or more are, in 
general, less useful than filters with lower ripple. Thus, 
it would appear that this awkward situation only occurs in 
filters which are not particularly useful. 

It should be noted, at this point, that the ripple 
figures for the high pass filter have been omitted. The 
reason for this is that the ripple contours are in exactly 
the same place as in the low pass filter. Of course, the 
definition of ripple, for high pass filters, is modified to 



r 



|H (f) I - lH(f„/2) 
i max ' I ' ' 



2 



(2-71) 



But the resulting contours are identical to the ones shown 
in Figure 2.57. 

The next figure (2.58) shows the iso-contours relating 
to the important quantity, "skirt slope". In analog filters, 
a well known rule of thumb holds that the skirt of the fil- 
ter (its transition region between passband and stopband) 
falls 6 db per octave for each real pole. Thus, this 
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Figure 2.58. Isocontours for the slope of the transfer 
function at cutoff for a second order low pass filter, 
a =1, a-|=2, d 2 =l. The h-plane triangle of stability. 
Slopes are given in terms of db per decade. 
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specification is often the determining factor for the order 
of Butterworth, Chebyshev, or elliptic filter. As can be 
seen in the figure / no such easy rule of thumb is possible 
in sampled analog filters. The slopes of the transfer 
function at cutoff are, of course, all negative (for low 
pass filters) . The greatest slope (absolute value) occurs 
for filters whose (b 2 fb^) pair is located in the lower right 
part of the triangle. This is also the region in which 
(Figure 2.55) cutoff frequency is extremely low. That is, 
these filters have a very narrow passband and extremely 
large slopes. The smallest slope contour is for -5. This 
contour cuts through many cutoff frequencies, from approxi- 
mately 0.23 f to the highest cutoff frequencies. Filters 

o 

of this nature have lower Q and are generally not so useful 
as the higher Q filters. If ripple is considered, we observe 
that for a maxim\im of 5% ripple, we have a wide range of 
possible slopes to choose from, or, we have a large range 
of cutoff frequencies from which to select. However, select- 
ing a cutoff frequency severely limits the possible skirt 
slopes . As can be seen by a brief study of the first three 
figures in this section (2.55, 57, and 58), a selection of, 

say, f = 0.15f and less than 5% ripple limits us to slopes 
c s 

between -51 and -38 db/decade. 

As an aid to locating the frequency for which the trans- 
fer function is a maximum. Figure 2.59 is included. The 
curves are once again equi-contours such that they are the 



209 




Figure 2,59. Isocontours for the position of the 
maximum of the transfer function for a lew pa s s second 
order filter, a =1, a =2, a =1. The b-plane triangle 
of stability. ° ^ 
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locus of points for which the maximum of the transfer 

function occurs at a given frequency. Naturally, if the 

ripple is zero, then the maximum of the transfer function 

occurs at f = 0. Thus, the zero ripple contour is also 

the contour for which f = 0. Since all the contours are 

max 

restricted to the region of non-zero ripple, the information 
on the figure is not so significant unless a designer is 
forced to use ripple. 

AN EXAMPLE OF DESIGN USING THE CHARTS (Figures 2.55, 
57,58,59) 

What is the frequency characteristic of the following 
transfer function in z? 

TT/\ 1 + 2Z^+Z^ 

H(z) = ^ s- (2-72) 

1 - 0.8z ^ + 0.5z 

Finding the frequency characteristic can be done graphically 
using the distances from zeros and poles, or can be done 
analytically by making the substitution 

z = (2-73) 

and solving the equation for various values of f. 

It can also be done, to a first approximation, by 
referring to the four previous figures. The following 
procedure is recommended; 

A. Locate the position of H_^^(f) from Figure 2.59 

max 

... In this case, f = 0.13f_. Place a point on a graph 

max s 
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with the coordinates (.13f , 1) which normalizes the 
maximum value of the transfer function to one. 

B. Consult the ripple chart (Figure 2.57) . 
Determine the ripple. In this case, r = 15%. Recalling 
the equation for ripple (Equation (2-68)), we then note 
that H(0) must be 0.7. Place this value on the graph. 

C. Since there is a minimum (zero) at f /2, 

s 

place a third point at the coordinates (0.5f , 0) . 

s 

D. Now locate the cutoff frequency from Figure 
2.55. In this case, f = 0.20f . Recalling the definiiton 
of cutoff frequency (Equation (2-69), we have H(f ) = 0.6. 
Place the point (0.2f , 0.6) on the graph. 

o 

E. Refer to the slope chart (Figure 2.58) . The 
slope at cutoff frequency is about -57 db/decade. This 

is fairly large. Draw a tangent line through the transfer 
function at cutoff frequency with approximately a -57 
db/decade slope. 

F. Finally, connect the points with a smooth 
curve. There will never be unexpected ripples in the fre- 
quency characteristic. The frequency characteristic will 
be monotonically increasing from dc to the frequency for 
which H is a maximum, and will then be monotonically 
decreasing from H to H(f /2). 

Figures 2.60 through 2.63 show the procedures for 
determining the approximate shape of the frequency transfer 
function. 
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points 



A SECOND EXAMPLE; 



Given a required specification (skirt slope, cutoff 
frequency, and ripple) , determine the appropriate values 
for and b 2 . 

The specifications to be used are 

= 0*2f 
c s 

ripple < 5% 

skirt slope maximum, but greater 

than 60 db/decade 

The following procedure is recommended; 

A. Locate the locus of points for which f = 0.2f . 
It is shown on Figure 2.55. 

B. Intersect this locus with the locus of points 
for which the ripple = 5% . All points on the cutoff fre- 
quency locus to the left of the intersection meet the first 
two specifications. 

C. Consult Fig\ire 2.58, the skirt slope chart. 

The -60 db/decade contour is not attainable with such a 
small ripple. The best that can be done is about -32 db/ 
decade. This occurs at b^^ = -0.62 and b 2 = 0.36. By 
actually going through this example, one can see that the 
skirt slope specification can be met, but only by changing 
one of the other two specifications. For lower cutoff fre- 
quencies, say f = O.lf , slopes of -40 db/decade can easily 

c s 

be attained at no higher ripple. 
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The use of recursive sampled analog filters using CCD's 
has been limited because of the stability problem and the 
anticipation that the parameters of the CCD may not be 
sufficiently constant to retain the desired filter specifi- 
cations. With this in mind, a sensitivity study was under- 
taken to determine how cutoff frequency and slope changes 
with marginal changes in CCD gain. From the canonical form 
of the second order filter, it can be seen that CCD voltage 
gain bears an important relationship to the coefficients a^^, 
b 2 . Usually these coefficients are implemented 
with potentiometers, that is, with passive devices whose 
maximum gain is 1.0. If we desire b^^ = 0.8, for instance, 
we must consider the gain of the CCD in cascade with the 
potentiometer for b^. The equation for b^^ is 

bjL = (K) (P) 

where K is the gain of the CCD and P is the potentiometer 
setting (see Section 2.2.1). If the gain of the CCD were 
1.5, then, for this case, we would have 

P = 0.8/1. 5 = 0.533 

It follows that a 5% change in CCD gain produces a 5% 
change in b^^. 
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For the second order coefficients, the gains of each 
CCD must be considered. The equations for these coefficients 
are 

^2 ^ ^2 ^ (2-74) 

^2 = K 2 P (2-75) 

Once again, it follows that an x% change in CCD gain produces 
an x% change in the effective coefficients. 

Figures 2.64 and 2.65 show the marginal change in cutoff 
frequency with respect to 5% changes in and b 2 (each 
separately). For Figure 2.64, the actual equation used 
for drawing the contours was 

f (1.025b,) - f (0.975b,) 

f = -£ f E (2-76) 

fc(b2> 

We observe in this chart that for most of the triangle, a 
change of 5% in b 2 causes a change of less than 2% in cutoff 
frequency, which should be tolerable for most applications. 
The particularly high change regions are near the triangle 
boundary for negative b^ and negative b 2 . Here the change 
in cutoff frequency is sometimes greater than 50%. However, 
a quick look back to Figure 2.55 where the cutoff frequency 
contours are plotted indicates to us that this is a large 
change in an extremely small number. Specifically, the 
region of maximum shift in cutoff frequency with respect 
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to a 5% change in b 2 is the region of smallest cutoff fre- 
quency (f < 0.05f ) . Thus, the absolute shift is probably 
c s 

small, but the relative shift is quite large. This region 
is probably only of academic interest anyway, since it is 
so close to the boundary of the triangle. Points close 
to the boundary are only marginally stable and should be 
avoided in practice. 

Figure 2.65 shows the same shift in cutoff frequency as 
a function of a 5% change in b^^. As may be seen from an 
examination of the figure, almost the whole inner region of 
the triangle above bj^ = -0.4 causes a change in cutoff 
frequency of less than 2%. Once again, the only region 
where extremely large variations are foxand is a region of 
extremely small cutoff frequency. 

It is understood that we have calculated the directional 
derivatives of a function f^(b 2 /b^) only in the directions 
of b 2 and b^^. That is, we have not calculated the contours 
for any other direction of variable change. Intuitively, 
we could determine the direction of maximiam change at each 
point in the triangle of stability and utilize it for design 
purposes. However, we believe that such a computation would 
only be of marginally more use than the present set of 
curves, and the present curves are easy to visualize for 
the design engineer. 

Example 2 required a given cutoff frequency, ripple, 
and skirt slope. If we were to include in the specification 
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that the cutoff frequency could not shift more than, say, 

5% with a 5% change in any coefficient, then we must resume 

our search for the proper coefficients. For the case where 

f „ = 0 . If „ , we find that the additional two charts offer 
c s 

new constraints on where the coefficients can be located. 

One possible point is b^^ = -0.9 and b 2 = 0.5. This point 
meets every constraint except the skirt slope constraint. 

The slope is only -57 db /decade. If this is not acceptable, 
and if no other specification can be changed in any way, 
then either the feedforward coefficients must be changed, 
or higher order filters must be designed. If the feed 
forward coefficients are changed, then the zero at f_/2 
may be moved or eliminated. For some applications, this 
is unacceptable. 

Finally, two other graphs have been plotted (Figures 
2.66, 2.67), each showing the shift in skirt slope with 
respect to a 5% change in the feedback coefficients. Pre- 
cisely the same procedure is followed to determine whether 
a region exists that meets all the specifications of the 
problem. Briefly, then, the overall procedure is as follows 

A. Using an overlay, trace the admissible regions 
for each constraint. 

B. Observe whether or not there exists any 
region in the triangle of stability which meets all con- 
straints. If so, select the appropriate coefficients and 
test. 
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shew percentage change in slope. ILov/ pass, second 
order filter, a =1, a, =2, a„=l. b-plane triangle 
of stability. ° ^ 
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Figure 2.67. Sensitivity of the slope at cutoff as a 
function of a 5% change in bp. Isocontours show 
percentage change in slope. '^Low pass, second order 
filter, a = 1, a = 2, a = 1. b-plane triangle 
of stability. ^ ^ 
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C. If not, then each constraint must be inspected 

closely to determine whether relaxation is possible. If 

not, then one can move the feed forward coefficients (a's) , 

or step up to higher order filters. Moving the a's has 

a deleterious effect on the zero at f„/2, while additional 

s 

order filters require more hardware and more study to 
determine the appropriate coefficients. 

2.4.4 DIRECT SAMPLED ANALOG FILTER DESIGN USING COST 
FUNCTIONS AND COMPUTER ITERATIONS 

One common engineering method for determining the values 
of a multidimensional vector has been the use of cost or 
loss functions. Briefly, one determines the required ideal 
function. Then one determines how much the relative cost 
is for deviation of the actual function from the ideal. 

This is usually sxibj actively performed by engineers knowl- 
edgeable in the field. A computer program is set up which, 
given a starting point, proceeds to "step" along the paths 
of the multi-dimensional vector until a value of the cost 
function is reduced to an acceptable level, or until a 
minimum is attained. The result is a set of values for 
each component in the multi-dimensional vector and a value 
for the cost function. Commonly, different starting points 
produce different "convergences" and several starting points 
must be used in order to insure that a global (or absolute) 
minimum of the cost function is reached rather than a local 
minimum. This is especially true for large order problems 
(where the multi -dimensional vector has many components) . 
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The first to utilize this method for the design of digi- 
tal filters was Steiglitz in 1970 [16], His procedure en- 
tailed defining a cost function (usually squared error, but 
not necessarily so) and a desired frequency response from 
dc to half the sampling frequency (the Nyquist frequency) . 

A starting point is chosen and the computer is used to 
locate points for the singularities which produce transfer 
functions which are successively closer to the ideal. When 

the computer finds a minimum (subject to small rounding errors 

« 

it stops and prints out the coefficients, the poles and zero 
locations, and the squared error. 

Using this method, it is not necessary to consult s-plane 
theory and documentation for an appropriate filter model. 
Although certain analytical proofs may be lacking (such as 
the proof that the elliptic filter provides the shortest 
transition region for any filter of equivalent order) , the 
results are generally so good that no improvement is necessary 
A program was generated as a part of this research to 
perform the fxinctions detailed above. Modifications to the 
program have been made to provide for proper weighting re- 
garding the seriousness of deviations from the ideal as a 

« 

function of frequency. As an example, it is often desired 
that a perfect zero be present in the stopband. This entails 
setting the coefficients a^, a^, and a 2 to constant values 
and letting the computer search for the optimum b^ and b 2 • 

Any nximber of different combinations can thus be used, each 
providing different weighting functions and producing a 
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different filter. The starting point for the direct design 
of a digital or sampled analog filter then, must be a fre- 
quency specification weighted for importance. No other 
input specification is necessary except for the order of 
the filter. 

Extensive use of the algorithm of steepest descent 
(Fletcher-Powell algorithm) suggests that several runs be 
taken for each problem starting at a different point each 
time. The well .known saddle point phenomenon occurs at 
times and different starting points will sometimes converge 
to different values. A comparison of the squared error 
function for each convergence will reveal the true optimum 
location of the singularities. 

One other pitfall is that the Fletcher-Powell algorithm, 
of course, makes no arbitrary limit on the position of 
poles and zeros. But our filter must have all its poles 
inside the unit circle of the z-plane in order to assure 
stability. Optimum locations of the singularities where 
the poles are outside the unit circle must be discarded 
since the filter is unstable. There is a way to eliminate 
such pole locations however. If we find that z^ is a pole 
location and jz^] > 1, then we multiply the transfer function 
by 



z 



z 




225 










I 



I 

I 

( 






which has a magnitude of when z is on the unit circle. 

This means that the inversion of a pole with respect to the 
unit circle does not affect the shape of the magnitude 
characteristic. Since the program is written to find optimum 
gain as well, such an inversion does not affect the char- 
acteristic at all. The program has been written to check 
automatically to see whether the poles and zeros are out- 
side the unit circle. If they are, they are inverted and 
the process is begun again at the new locations of poles 
and zeros. 

In addition to being a replacement for the more familiar 
design by analog models, this process works for any arbitrary 
frequency characteristic specification. Thus, we are not 
limited to low pass and high pass models. It is perhaps 
interesting to note that our frequency characteristic is 
defined for discrete frequencies of interest. They can, 
of course, be as close together as desired, but since the 
computer is digital, the characteristic cannot be totally 
free of the possibility of unwanted excursions BETWEEN 
specification points. However, as long as the niimber of 
specifications greatly exceeds the order of the filter, this 
will not be a problem. On the other hand, it can be shown 
that if the order of the filter is equal to twice the number 
of frequency specifications, then the specification will, 
in theory, be met exactly. (This is a result of having 2n 
degrees of freedom for nth order filters.) However, to 
meet each frequency characteristic specification precisely 
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runs the risk of having excessive deviation between the 
frequency specifications. Actually, this problem is more 
academic than practical, since in the overwhelming majority 
of cases, the order is low and the frequency specification 
is given as a continuous function. It is well to remember 
that, in the case of a piecewise constant transfer function, 
more than one frequency, in general, should be included as 
a specification for each constant region. 

Figures 2.68, 2.69, and 2.70 show some of the results 
of using the Fletcher-Powell-Steiglitz method with the 
program generated during this research. Figure 2.68 shows 
a second order low pass filter where the feed forward 
coefficients were permitted to change. Figure 2.69 shows 
the same type filter where the feed forward coefficients 
were held constant at (1,2,1). Figure 2.70 shows the first 
"tooth" of the filter in Figure 2.69. 

Section 2.5 Experimental Results and Theoretical Calculations 
of Representative Second Order Comb Filters 

During the design of the second order comb filter using 
two CCD's we dealt with several complexities not found in 
the first order case: 

A. Two CCD's had to be found which had similar 
characteristics in order to "share" bias levels. One of 
our requirements was that these bias levels be common for 
corresponding gates on each CCD. We tested and characterized 
several CCD's in the search for such similarity. We finally 
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Figure 2.68. Transfer function of a second order, low pass.sampled analog filter where 
feedforward coefficients were permitted to change. 
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Figure 2.69. Transfer function of a fourth order* low pass, sampled analog filter 
where feedforward coefficients were held constant at (1,2,1). Notice the zeros at 
multiples of half the sampling frequency. 
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Figure 2.70. First "tooth" of the transfer function of a second order, low pass, sampled 
analog filter whose feed forward coefficients were held constant at (1,2,1). Notice 

the zero at half the sampling frequency (010 on the horizontal axis). 



decided on two CCD's from the 747 lot, made by TRW. The 
operating curves which seemed best were curves for the input 
on gate 3 and making serve as the constant parametric 
voltage. During the course of the experiments, we constantly 
monitored the actual output waveform of the CCD since they 
tended to drift off the proper bias levels, and also for 
filtering reasons since feedback often causes large signals 
to appear internally. 

B. More hardware was needed. Instead of needing 
only one sample-and-hold device, we initially used two . . . 
one to restore the analog signal after the first CCD, and 
one to restore the analog signal after the second CCD. This 
did not turn out to be adequate, and another sample-and-hold 
device was inserted at the input in order to be sure that 
theoretical cancellation would take place. That is, at the 
summers, we insured that all THREE input signals were of a 
sampled-and-held form. Also required were additional level 
shifters. Naturally, the output level had to be shifted so 
that the ac waveform out had an average value of zero. This 
was done by using two "up" shifters following each CCD. 

Also necessary were two "down" level shifters since the 
input of each CCD had to be shifted substantially further 
negative in order for the proper bias conditions to exist. 
These level shifters were built from operational amplifiers 
with high input impedance. As a matter of fact, we had 
improved the input impedance to the summers and level 
shifters so that we did not have to make the "correction 
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factor" discussed in Section 2.2.1. The only correction we 
had to make in the setting of the potentiometer coefficients 
was for the voltage gain of the CCD. 

Figure 2.71 shows the physical arrangement of circuit 
components for the evaluation of the second order sampled 
analog recursive comb filter. 

We planned the experimental work by investigating Figures 
2 . 55 through 2 . 59 for several transfer functions which would 
show feasibility of sampled analog recursive filters and 
also illustrate the power of the b-plane triangle of stability 
concept. Table 2.5 shows our complete experimental plan. 

Some points on the b-plane were selected for their monotonic 
transfer functions, while others were chosen for their 
large ripple. All the filters tested were low pass, but 
since the high pass triangle contours are precisely the 
"upside-down, mirror image" of the low pass triangle con- 
tours, we believe that little generality was lost. Some 
pairs of b’s were selected which were close to the boundary 
of the triangle, and their transfer functions exhibited, 
in general, higher selectivity (Q) , than those chosen 
nearer the center of the triangle (the origin) . 

There is one difficulty in setting coefficient values 
which should be passed on to those who will make further 
measurements. For that reason, it is included here: 

A. Settings of b^ may range from -2 to +2. One 
good way of deriving the proper potentiometer value is to 
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TABLE 2.5 Data Taking Plane 
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t calculate the recursive equations for each output in time as 
a function of an input unit pulse or unit step. Then, 
while observing the output waveform on the oscilloscope, 

, impress either the unit pulse or the unit step at the input. 

I The unit step input was found to be more convenient. Since 
a first order filter is simpler, b^ may be set while a .2 and 
I b 2 are set at zero. This procedure was followed until the 
magnitude of b^^ exceeded unity. At this point, a first 
order filter can no longer be designed with such feedback, 

I and the other feedback coefficient, b 2 must be introduced 

I 

! into the design. 

The recursive equations we used for setting coefficients 
whose magnitude exceeded unity are shown in Table 2.6. The 
experimental procedure is to impress the unit step at the 
input, observe the output, and then change the potentiometers 
until the proper output is reached. Record the settings, 
and proceed with the experimental work. 

Problems were found in the cases where the pair (b 2 /b^) 
were close to the boundary of the triangle. At times, no 
data at all could be taken since the filter was unstable. 

I Four cases where no data at all was taken are shown in 
Table 2.5. 

Figures 2.72 through 2.89 are tabulated here as represen- 
tative of the data taken during the course of this experiment. 
All the graphs are different from their coi^terparts in the 
first order case because of the change in the y-scale. It 
is logarithmic. The "wiggles" in many of the curves occur 
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Figure 2.72. Theoretical and experimental transfer function for a second order 
recursive cont> filter, b = 0, = -0.7. Dotted lines are experimental points. 
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Figure 2.73. Theoretical and experimental transfer function for a second order 
recursive comb filter, b = 0, b^ = -0.3 Dotted lines are experimental points. 
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Figure 2.74. Theoretical and experir.iental transfer function for a second order recursive 
comb filter, b, = 0, b = 0.3. Dotted lines are experimental points. 
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Figure 2.75. Theoretical and experimental transfer function for a second order recursive 
comb filter, b, = 0, b = 0.7. Dotted lines are experimental points. 
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Figure 2.76. Theoretical and experimental transfer function for a second order 
recursive comb filter, b, = -0.3, b = -0.3. Dotted lines are experimental points. 
(LOG SCALE) ' 2 



H(f) 



O 




ZC^-' 



zzrj- 






241 



Figure 2.77. Theoretical and expeririental transfer function for a second order 
recursive comb filter, b = -0.3, b- = +0.3. Dotted lines are experimental points. 
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Figure 2.78. Theoretical and experimental transfer function for a second order recursive 
comb filter, b = -0.3, bo = +0.7. Dotted lines are experimental points. 
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Figure 2.79. Theoretical and experinental transfer function for a second order recursive 
comb filter, b, = +0.3, bo = -0.3. Dotted lines are experimental points. 
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Figure 2.80. Theoretical and experimental transfer function for a second order 

recursive comb filter, b = +0.3, b = +0.3. Dotted lines are experimental points, 
fino 1 2 
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Figure 2.81. Theoretical and experimental transfer function for a second order recursive 
comb filter, b = +0.3, b^ = +0.7. Dotted lines are experimental points. 
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Figure 2.82. Theoretical and experimental transfer function for a second order 
recursive comb filter, b = -0.6, b = +0.1. Dotted lines are experimental points. 
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Figure 2.83. Theoretical and experimental transfer function for a second order recursive 
comb filter, b, = -0.6, b = 0.7. Dotted lines are experimental points. 
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Figure 2.84. Theoretical and experimental transfer function for a second order 
recursive comb filter, b = +0.6, b- = +0.3. Dotted lines are experimental points. 
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Figure 2.85. Theoretical and experimental transfer function for a second order recursive 
comb filter, b, = -1.0, b^ = +0.3. Dotted lines are experimental points. 
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Figure 2.86, Theoretical and experimental transfer function for a second order recursive 
comb filter, bi = -1.0, b = +0.7. Dotted lines are experimental points. 
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Figure 2.87. Theoretical and Experimental transfer function for a second order recursive 
comb filter, b, = +1.0, b„ = +0.7. Dotted lines are experimental points. 
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Figure 2.88. Theoretical and experimental transfer function for a second order recursive 
comb filter, b-, = -1.3, b^ = 0.45. Dotted lines are experimental points. 
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Figure 2.89. Theoretical and experimental transfer function for a second order 
recursive coni) filter, b, = +1.3, b = +0.8. Dotted lines are experimental points. 
(LOG SCALE) ^ 
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used to set bj^ when b^^ = 1. 
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' 0, b^ =1, bj^ = 1, i:>2 — 1 


TIME 


INPUT 
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(l+bj^)bj^+l+b2 = l+b^+bj^^+b2 
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1 
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at a low response level and are the results of the system 
noise as well as imperfections in the filter itself. This 
was not apparent in the first order experiments because of 
the linear scale. In many cases agreement between theory 
and experiment was extremely good . . . notwithstanding the 
somewhat erratic results seen at low output level and at 
frequencies near the Nyquist limit (10 kHz) . 

All the theoretical curves are corrected for the sample- 
and-hold effect and also for the inefficiency of the device 
which we conservatively estimated at 0.001. Thus, the 
theoretical curves are somewhat different from what one might 
expect. The inefficiency causes the null (at f /16, 3f /16, 
5f /16, ...) to be progressively poorer as frequency is 
increases. It also causes the peaks of the output to change 
with frequency and sometimes , we get a triple peak as in 
Figure 2.88. The triple peak would not be possible if 
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inefficiency were zero since, in that case, the values of 
the feedforward coefficients would always be precisely 
(1,2,1) so that both zeros in the transfer function would 
lie at z = -1 on the z-plane. But the effect of ineffi- 
ciency moves the zeros away from that point, thus causing 
multiple peaks in the passband and a deterioration of the 
null in the stopband. 

Of course, the sample-and-hold effect is a familiar 
sin x/ X multiplicative effect. 

Once again, as we did in the experiments in first order, 
we find that several transfer functions have the disturbing 
tendency to rise with increasing frequency. We believe 
that the major factor contributing to this rise in the 
increase in reactance causing a decrease in b^^ and b 2 at 
higher frequencies. A quick look at the individual experi- 
ments which manifested this phenomenon confirms that the 
rising occurred mainly when the feedback was positive; 
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Since these same potentiometers were used in previous experi 
ments as well, it is appropriate to determine whether or not 



others have found this phenomenon. In a thesis completed 
here in June, 1975, Ensign B. Freund determined that the 
RETICON (SAD) serial analog delay device exhibited similar 
properties when used in first order filters. One common 
set of equipments used in the different experiments were 
the potentiometers. 

In order to investigate this increase further, we backed 
off from our assumption of an equivalent second order system 
and took the more conservative view that our device is a 
sixteenth order filter. Further, we used a set of coeffi- 
cients (b 2 fbj^) which, when 7.8 kHz was reached (at the 
input) , caused the filter to become unstable. Now this 
behavior cannot be predicted by the second order transla- 
tion, but it can be easily \anderstood by a z-plane plot of 
the sixteen poles, which is shown in Figure 2.90. 

In order to determine that this was the true situation, 
we advanced the frequency to about 9.9 kHz where once again 
the filter exhibited instability. This instability per- 
sisted until past 10.0 kHz and finally subsided at 10.15 
kHz. Increasing the frequency further (now beyond the 
Nyquist rate), we found the filter to be stable until 12.2 
kHz was reached. Here once again we observed unstable 
operation. Our sixteenth order theory told us that we 
should be prepared for an instability at a frequency which 
corresponded to the complex conjugate of the instability 
at 7.3 kHz. The frequency at the complex conjugate is 12.2 
kHz . 
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Figure 2.90. The pole plot for our sixteenth order 
filter which displays instabilities at 7.8 Khz and 
9.9 Khz. 
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Going further, we increased the frequency to 20 kHz, 
but no other unstable conditions appeared. It was postu- 
lated that none would appear until we had traveled all the 
way around the unit circle (on the z-plane) to the frequency 
corresponding to 7.8 kHz + one revolution or 27.8 kHz. We 
found this unstable condition at 28.0 kHz. Thus, we found 
that our filter, although not conforming to second order 
translation theory, did conform to sixteenth order filter 
theory. Of course, there was no way to measure the radial 
location of the pole at the point corresponding to 7.8 kHz, 
since the pole may lie outside the unit circle. 

SUMMARY 

Out of 33 planned experiments, we were able to complete 
29. The other four proved so unstable that no data could 
be taken. Of the 29 sets of data, we find that in many 
cases, close agreement between theory and experiment existed 
at least for the first two "teeth." Also, we found that 
the shapes of the transfer functions were very close (once 
again, for the first two teeth) to that predicted for them 
by inspecting the b-plane triangle of stability data. 
Measurements of cutoff frequency and ripple, especially, 
were compared to their predicted values with close agreement 
in most cases . These measurements were made on the first 
"tooth" of the transfer function which is, of course, the 
"best behaved" region in the frequency domain. No attempt 
was made to deal with the standard z or bilinear z transforms. 
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The standard z transform, once again, produces (for second 
order) a transfer function which has zeros only at infinity 
or zero. The bilinear z transform does produce, in general, 
two zeros and this was, in fact, what was used. However, 
rather than use the z-transform to design our filters, we 
obtained our coefficients directly from the b-plane triangle 
of stability. Figure 2.91 shows a three dimensional look 
at cutoff frequency as a function of (b 2 /b^) . 

Section 2.6 Conclusions 

Our goal when undertaking the experimental work was 
to show that sampled-analog devices, such as the CCD, were 
appropriate to the design of recursive comb filters. At 
the outset, we were one of the only groups working on this 
project, although many had utilized the charge transfer 
devices for delay elements in transversal filters. On the 
basis of the data contained in this chapter, we believe that 
we have shown that sampled analog devices can be used as 
the delay element in recursive comb filters. 

There are several other conclusions which are significant 

A. The standard z transform is not very useful 
when designing to a frequency specification. 

B. The bilinear z transform is much more useful, 
since it provides that the range of tlie sampled analog 
filter is precisely equal to the range of the analog filter 
from which it was derived. 
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ire!. 91. Cutoff frequency (vertical scale) as a function of (b^.b-i). The base 
ihe)-plane. The bo axis runs from foreground (b„ = -1 ) to background (b„ = 1). 
ti^ixis runs across the page from +2 to -2 (left^to right). ^ 
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C. One can make the eight-to-first order and 
sixteenth-to-second order translations for a specific set 

of feedback and feedfoirward coefficients, thus affording the 
engineer with a powerful tool with which to design filters. 

D. Inefficiency, the sample-and-hold effect, and 
the frequency response of the potentiometers (and the 
operational amplifiers, too) detract from perfect operation 
of the filters. 

E. Second order sampled analog filter design is 
best done with a general purpose computer, using the program 
this research generated, or its equivalent. 

F. Analog to Sampled analog second order filter 
design is not very tractable, and is not any better (in 
fact, worse) than using the computer program. 

G. For the engineer at the bench, one easy way 

to predict the changes in the transfer function as a result 
of modifications of the coefficients is the b-plane triangle 
of stability curves. Cutoff frequency, ripple, the posi- 
tion of H , and the slope at cutoff have been calculated 

max ^ 

and isocontours are now available. Also, marginal changes 
in these quantities (especially cutoff frequency and slope) 
as a function of small incremental changes in b^^ and b 2 are 
now established and shown in four of the triangle drawings. 

H. The eight bit CCD does not provide enough 
"teeth" at frequencies far below fg/2 to make a true comb 
filter. The RETICON SAD device is better. Still better, 
at least insofar as number of delay stages is concerned. 
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is the Fairchild CCD {455 delay stages) recently 
acquired. 

I. The magnitudes of the ''null" frequencies 
in first and second order filters are about the same. 
Therefore, if this were the dominant specification, one 
would be advised to use first order filters for simplicity. 
On the other hand, second order filters make it possible 
to achieve higher selectivity, Q. 
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CHAPTER THREE 



APPLICATION OF THE SAMPLED-ANALOG RECURSIVE COMB FILTER 
3.0 Introduction 

The contents of this chapter serve to supply the link 
between device/circuit research and systems applications. 

The primary application we have selected arises rather 
naturally in radar systems, that of the elimination of 
the echoes of stationary targets and the enhancement of the 
echoes of moving targets. Such a process is usually per- 
formed by a device called a Moving Target Indicator (MTI) . 

The use of MTI is not novel . . . navy air search radar systems 
have had designed-in MTI for twenty years. The problem was, 
however, that cancellation was generally done (not too well) 
by analog means and by use of a device called a quartz delay 
line. The physics of this device are adequately described 
in Skolnik's famous text. Introduction To Radar Systems . 

The following subjects will be addressed; 

Section 3.1 . Theory of echo cancellation; cancellers, 

filters 

Section 3.2 . Time domain truncation; effects on 
echo spectrum and cancellation 

Section 3.3 . The recursive comb filter as the MTI 
processing element 

Section 3.4. Conclusions. 
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Section 3.1 Theory of Echo Cancellation; Cancellers, Filters 



Echo cancellation and the distinguishing of stationary 
target echoes from moving target echoes is made theoretically 
possible by the physical phenomenon of doppler frequency 
shift. Briefly, doppler frequency shift occurs between a 
transmitted wave and a received echo wave when there is 
relative motion between transmitter and target. A discussion 
of doppler frequency shift, which can be observed not only 
in electromagnetic waves, but also in acoustic waves, can 
be found in most basic physics textbooks. If the relative 
motion between the two objects serves to decrease the distance 
between them, then the doppler frequency is considered posi- 
tive, and the frequency of the received signal is higher than 
the transmitted frequency. Conversely, if the relative 
motion between the two bodies serves to increase the dis- 
tance between them, then the doppler frequency is negative 
and the frequency of the received signal is lower than the 
transmitted frequency. If the distance between the objects 
is unchanged during the time of the transmission and recep- 
tion of the pulse, then the frequency of the received signal 
is exactly the same as the transmitted frequency and the 
doppler frequency is said to be zero. 

Since echoes from stationary and moving targets occupy 
different regions in the frequency domain, it is theoretically 
possible to distinguish between them ... to build equipments 
which will ignore stationary targets, but sound alarms on 
fast moving targets. After the echoes have been categorized 
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into stationary and moving targets, of course, the system 
can be made to produce almost any response . . . including the 
automatic firing of a missile. But it is the initial problem 
of determining the radial speed of a target which has proven 
difficult over the years. 

The doppler frequency is given by the following equation: 



f 



d 



^ ’'r 



(3-1) 



where v^ is the radial speed of the target, f^ is the trans- 
mitted frequency, and, of course, c is the speed of light. 

It perhaps provides insight here to determine real doppler 
frequency for real parameters. The transmitted frequency 
of the AN/SPS-29 series air search radars is about 215 mHz. 

A target traveling at 400 knots of radial speed will produce 
a doppler frequency shift in this radar set of 296 Hz. 
Observe that the doppler frequency is about 0.000138% of 
the transmitted frequency. Since such a small percentage 
could rarely be detected by RF discriminators, the usual 
procedure is to shift the entire received echo spectrum to 
baseband, or so that the transmitted frequency is shifted 
to zero frequency. Then a bandpass filter could be used to 
distinguish the target echo at 296 Hz. 

The analysis of doppler frequency has, of course, been 
done . Skolnik proceeds as follows : 

A radar pulse is transmitted for t < t^ + T and 

is given during that time by 
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(3-2) 



= A sin 2-irf^t 



The received signal is 



''echo = ® sin(2n(£^+f^) 



4irf 



(3-3) 



In the argvunent of the first term represents the 

shift in frequency due to the doppler effect. The second 
term represents the shift, relative to in phase of 

the echo due to the target's range. 

Maximum unambiguous range is R = cT/2 where T is the 

1113. X 

time between pulses . 

We extract the doppler information by mixing (multiplying) 
the t)y ths ^echo* want the difference frequency 

which will be of the form 



V = C sin (2TTf ,t - d>) 
o a ^ 



(3-4) 



Now there have been historically two major ways of processing 
video, or of presenting echo information to the operator. 

The first, chronologically, was the A-scope which provided 
to the operator a pulse at some distance along the horizon- 
tal axis which corresponded to the range of the target. MTI 
processing for the A-scope took the form of the well known 
"Butterfly" effect. Consider Equation (3-4). If the radial 
speed of the target is zero, then the doppler frequency is 
zero. Then we are left with 
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(3-5) 



= C sin {-(j>) 

which is a constant for all time. 

On the other hand, if the radial speed was not constant, 
then f^ 5 ^ 0 and oscillates in time. Thus, the operator, 
on observing an oscillating target echo on his A-scope knew 
with no further processing that the target was a moving 
target. On the other hand, he could safely ignore constant 
sized pulses since they represented stationary targets. 

Figures 3.1 and 3.2 show, respectively, the simulated A- 
scope presentation of a stationary target echo and a moving 
target echo. 

The second form of display has been the popular PPI 

(plan position indicator) which is shown in Figure 3.3, 

and is well known to almost every naval officer. Here the 

plot is in terms of angle and range, not amplitude. Thus, 

the operator cannot distinguish between a c onstant- amplitud e 

pulse and an oscillating-amplitude pulse, since the indication 

on the PPI is in terms of power. In order to provide some 

stationary vs. moving target processing, it is necessary to 

pre-filter the echoes prior to introduction to the PPI. 

The simplest method is the one zero canceller shown in Figure 

3.4. Here the first echo from a target arrives at time t^, 

and the second at time t + T. These signals are introduced 

o 

into the subtractor with the following results: 
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Figure 3.1. A-scope > 
presentation of a stationary 
target. Notice the constant 
amplitude. 




Figure 3.2. A-scope presentation 
of a moving target. Notice the 
so-called "Butterfly" effect. 




Figure 3.3. The plan position indicator 
(PPI) display. 
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= A sin (2TTf^t - ())) 



V ^2 = A sin (2Trf^(t+T) - (fi) 

^ (2Trf^(t+j) - 4>) (3-6) 

As can be seen by an inspection of Equation (3-6) , if the 
doppler frequency is zero, then the amplitude of the response 
for all time is zero. On the other hand, if the doppler 
frequency is not zero, then there is some nonzero amplitude 
response which is passed to the PPI. 

It must, of course, be pointed out that there are many 
different doppler frequencies for which Equation (3-6) is 
zero for all time. Such doppler frequencies are given by 

f^ = n/T ; n = 0,1,2,... (3-7) 

The meaning of this last equation is that if the target has 
a radial speed such that the doppler frequency satisfies 
Equation (3-7) , then this processor will provide zero response. 
The physical interpretation is that the target will never 
be detected. Such speeds are called BLIND SPEEDS and many 
different designs have been tried to eliminate them, or at 
least reduce their effects. To determine how significant 
these blind speeds may be, an example with real parameters 
is in order. Referring to the aforementioned AN/SPS-29 
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iigg', 



series air search, we postulate the following characteristics: 



= 215 mHz 

PRF = 1/T = 266 Hz 

R = 250 mi. 
max 

For these parameters, the first blind speed occurs at 

= c/2f^T = 305 knots 

The first blind speed, then, is 305 knots, a speed well 
within the operating capability of many aircraft and missiles. 
As a matter of fact, the second and third blind speeds, 

610 and 915 knots, also must be considered as possible target 
speeds. The significance, then, of blind speeds is great. 

It must be emphasized that no amount of pulse shaping 
can eliminate the undesirable zero response at blind speeds. 
This phenomenon is caused by the periodicity of the target 
echoes, which, in turn, is caused by the periodicity of the 
radar transmitted pulse. One way of looking at the phenomenon 
which is appropriate to our discussion, is that the radar 
echoes represent a sampling of a signal waveform of the 
doppler frequency. Figure 3.5 is drawn to make this aspect 
of the phenomenon clearer. 
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Figure 3.5. Radar encounters n target whose radial speed is such that 
the doppler frequency is exactly equal to the PRF. Since all target 
echo amplitudes are equal, the target provides the same input pulse train 
to the MTI as a stationary target. The response of the MTI filter is zero 
at this BLIND radial speed. 




Figure 3.6. Radar encounters the same target (equal doppler frequencies) 
but with staggered PRF, the target echo amplitudes are not equal (as shown). 
Thus cancellation will not take place and the response of the MTI filter 
will be nonzero 
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One way of eliminating blind speeds is to change the 
pulse repetition frequency. Several methods of change have 
been proposed and used. One is called "staggered PRF" which 
utilizes two different PRF's and shifts them after each 
series of pulses, or sometimes, after each pulse. Complex 
gating in the radar set insures that the targets at equal 
ranges will be processed together in the MTI. The other 
method is called "Jittered PRF" which selects, at random, 
from a range of PRF's, a new PRF for each pulse. This method 
also requires complex gating to insure comparison of target 
echoes at appropriate ranges. The effect of staggered PRF 
is shown in Figure 3.6. The doppler frequency is drawn for 
all time and is sampled at those times during which it is 
illuminated by a radar pulse. Notice that the amplitudes 
of the target (traveling at a blind speed) echo are not 
equal, and thus, will not cancel completely in MTI canceller. 

As a matter of fact, the first blind speed now occurs at 
the first multiple such that 

n/T^ = m/T 2 where n and m are positive integers. 

T^ and T2 are the pulse repetition Intervals respectively. / 

The traditional means of accomplishing analog delay has 
been the acoustic delay line. The steps involved in using 
this delay line are 

1. Convert electromagnetic energy to acoustic energy. 
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2. Insert the acoustic energy into a quartz or 
mercury delay line . . . which is a crystal shaped so that 
the energy will bounce from one face to another in a pre- 
determined manner in order that it will appear at the output 
after a specified delay time. 

3. Reconvert the acoustic energy output to electro- 
magnetic energy. 

The disadvantages of the analog delay have been known 
for sometime. They are 

1. The delay line and associated transducers are 
quite lossy providing severe insertion loss . 

2. The processor is bulky and large. 

3. The crystal is sensitive to temperature varia- 
tions . Changes in temperature change the shape of the crys- 
tal slightly causing a shift in the nominal delay time. 

Such shifts in delay time destroy proper operation of the 
canceller. 

4 . There are bandwidth limitations for some 
applications . 

5. The crystal is sensitive to shipboard shock 
and vibration. 

6. It is difficult to manufacture the crystal to 
required specifications ... or if it is done, it is extremely 
costly. Slight imperfections can cause secondary responses 
(ghosting) . 

7. Inhomogeneities in the material itself can 
cause spurious responses as well. 
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Today, many radar systems have digital processors. 

The raw analog target data is first converted to digital 
data for insertion into digital processors. Digital filters 
have been built and are utilized in MTI applications. Even 
some recursive filters are utilized in this way. Of course, 
in any analog to digital system, there are certain disadvan- 
tages ; 

1. Analog to digital conversion produces some 
inevitable quantization noise ... so-called noise resulting 
from the rounding off of the actual analog values. This 
quantization noise takes on considerable significance when 
using recursive filters since, in recursive filters, the 
same digital data is fed back to the input repeatedly thus 
increasing the quantization error. A brief example will 
show what we mean. Suppose we have a feedback coefficient 
of +0.9 and we insert a unit pulse into the filter: Round 

off the output to the nearest tenth. 



Start with input at time 


zero = 


1.0 


. . . output = 1.0 


time 


one = 


0.9 


. . . output = 0.9 


time 


two = 


0.81 


. . . output = 0.8 


time 


three= 


0.72 


. . . output = 0.7 


time 


four = 


0.63 


. . . output = 0.6 


time 


five = 


0.54 


. . . output = 0.5 


time 


six = 


0.45 


. . . output = 0.5 


but no matter how much further we go 


, we 1 


cannot change the 



output from 0.5. In this extremely simple example, we have 
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shown that the quantization error can be quite significant 
in digital systems. The steady state error in our system 
was the difference between zero and half the original value 
of the signal, clearly unacceptable. 

Of course, ■ there are known ways of combatting such catas- 
trophic errors. The first is to expand the number of quan- 
tization levels. As more quantization levels are provided, 
the error is consequently reduced. Of course, this means 
more hardware, and greater complexity of design. Perhaps 
of greater import than the actual hardware is the computer 
study which must be undertaken to determine exactly how 
many more quantization levels are required. Several theses 
and papers originating at the Naval Postgraduate School 
have led the field in this subject (see especially, work 
by Dr. S. R. Parker) . These computer studies must be made 
by radar/processing engineers and, thus, are costly. 

Another way of reducing the effect of digital quantiza- 
tion error is to prevent the samples from being multiplied 
over and over by the same quantized gain coefficient. This 
takes the form of eliminating the canonical for digital 
filters and substituting transversal filters, or, at least, 
using cascaded recursive sections of not more than two 
orders each. Quantization error in digital systems can be 
reduced, but it cannot be eliminated, and we are just 
beginning to understand methods to predict significant 
quantization error and take countermeasures. 
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2. Making the conversion from digital back to 
analog is also a source of noise. 

3. Naturally, the hardware involved in each of 
these conversions increases the complexity of the equipment, 
and also increases the design time of the radar system 
itself . 

There is one other important factor to notice about 
these cancellers. Since they are filters, they have a 
transient response to any input signal. As with all linear 
filters, their frequency transfer functions are the Fourier 
transforms of their unit pulse responses. A more useful 
specification for these cancellers is the so-called indicial 
function which is nothing more than the response of a linear 
system to a unit step input. The relation between the unit 
pulse response and the indicial function is 

h(nT) = A(nT) - A((n-1)T) (3-8) 

where h is the unit pulse response and A is the indicial 
function. Ideally, then, we would want the indicial func- 
tion, A, to be identically zero since the echo from the sta- 
tionary target is tantamount to receiving a unit step input. 
Unfortunately, the indicial function is not zero for all 
time, but has a nonzero transient. This transient, of course, 
can be calculated for any given filter, but it is often 
erroneously ignored, and it must be considered when assessing 
the value of a filter. We have measured some of the time 
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responses in the laboratory of various cancellers. The 
measurements agree quite closely with the predicted values 
of h and A. Figure 3.7 shows two cases of the use of filters 
which have zeros in their transfer function at dc, and yet 
the indicial function and unit pulse response are nonzero 
for a substantial time period. For the case in Figure 3.7a, 
the feedback coefficient of our first order filter was only 
-0.3 and the indicial function dies out rather rapidly. The 
magnitude of the function falls below O.IA(O) at n = 5 and 
stays below that value. On the other hand, the filter with 
the feedback of -0.9 has an extended response such that A 
does not fall below O.lA(O) until the 29th pulse output. 

Naturally recursive filters are not unique in having 
nonzero transient behavior. Transversal filters (cancellers) 
also exhibit this troublesome characteristic. In the can- 
celler shown in Figure 3.4, the transient response to a unit 
step input can easily be calculated and is shown in Figure 
3.8. The transient response in this case is minimal, but, 
of course, it is accomplished at the cost of a nonideal 
frequency characteristic which is shown in Figure 3.9. The 
passbands of the filter are clearly not flat, a circiamstance 
which can, at times, lead to a loss of detection. This can 
arise when the operator notices that the level of the video 
is too high for a target traveling at a radial speed for 
which the response of the filter is maximxam. The operator's 
impulse is to reduce the gain of the radar set, and when he 
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Figure 3,7a» Unit pulse response of first 
order recursive filter, = 0.3 




Figure 3.7b. Unit pulse response of first order 
recursive filter, b^ - 0.9 
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input = unit step 



t ti me 

0 



indicia! function = response to unit step input 



time 



Figure 3.8. Transient response of one zero-no pole MTI canceller, 
to unit step input. Filter v/ill exhibit another pulse when 
input falls to zero. Filter is shown in figure 3.4. 
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Figure 3.9. Transfer function of first order (one zero, no pole) canceller for MTI. 



does so, he eliminates the possibility of detecting targets 
at some of the marginal radial speeds. Flatness in the pass- 
band can, of course, be better achieved by increasing the 
order of the filter, but this implies lengthening the unit 
step response. As can be seen from Chapter two, and par- 
ticularly the second order case, filters close to the ideal 
frequency characteristic can be designed. Filters whose 
time characteristic is close to the ideal can also be de- 
signed. But filters "close” to both, cannot. 

t 

Section 3.2 Time Domain Truncation: Effects on 

Echo Spectrum and Cancellation 

There is one deviation from theory which is so signi- 
ficant in radar systems that it deserves a section all by 
itself. It is not often considered in designing a radar 
system with MTI, but its effects are far reaching and 
important. It is the time domain truncation that almost 
every search radar exhibits . 

Consider a normal air search radar. It is mounted high 
on a ship and rotates in azimuth. Pulses are transmitted 
and are returned from targets before the rotation of the 
antenna can move far enough to ^avoid receiving the echoes. 

If this were not so, no echoes would ever be received. The 
significant development is that the antenna only stays "on 
target” for a small finite time. It will be shown that 
this small finite time is extremely important in determining 
the response of the MTI and, thus, the ability of the radar 
to distinguish stationary from moving targets. 
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The small finite time of "on target" enables the radar 
to send out k pulses, each of similar width. It is this 
number, k, which turns out to have such a large effect 
on the spectrum of the target echo train and the subsequent 
cancellation. 

As k approaches infinity, of course, we have the ideal 
case which we had assumed in Section 3.1. The spectrum of 
the stationary target is periodic (spikes at dc, and each 
multiple of the PRF) , . The spectrum of moving targets is 
also periodic with period equal to the PRF. Clearly, the 
filter with the frequency characteristic having zeros at 
dc and the multiples of the PRF will eliminate the stationary 
target completely and pass the moving target echoes. Figures 
3.10a and 3.10b show the ideal frequency spectrxim for 
stationary and moving targets. 

As k is reduced from infinity, we no longer have the 
ideal case. Physically, this represents the antenna scan 
increasing or the pulse repetition frequency decreasing. 

The spectrum of the input, then, can no longer be discrete, 
although it will still be periodic. This follows from the 
familiar properties of the Fourier transform pairs; If 
a function is not periodic in one domain, then it will be 
continuous in the other domain. 

Since the spectrxam has become continuous (and if it is 
not everywhere zero) , then there must be a band of frequencies 
in which substantial signal energy lies. If we choose some 
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Figure 3.10a. Ideal spectrum of the target echoes from a 
stationary target. Spacing of the lines in the spectrum is 
1/T (PRF). 




Figure 3.10b. Ideal spectrum of the target echoes from a 
moving target. Target radial speed is such that the doppler 
frequency is f^ = 1/3T. The spacing of the lines is 1/T (PRF). 
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k, we can calculate the frequency characteristic directly. 

We have selected three values of k to demonstrate the 
effect of time domain truncation: For k = 2, Figure 3.11 

shows that the energy in the signal (still from a stationary 
target) is hopelessly spread throughout the band from dc 
to 1/T. For k = 6, we can see in Figure 3.12 that the 
energy is confined to the major lobes at dc and 1/T, but 
still a s;abstantial amount of energy exists in the nominal 
passband of some MTI filter. For k = 10, Figure 3.13 shows 
that the level of the minor lobes has fallen from the pre- 
ceding case and that most of the energy lies in the region 
around dc and 1/T. The relative amount of energy of the 
signal in the nominal passband of the MTI filter compared 
to that in the stopband can, of course, be calculated for 
any given filter and any given number of "hits on target". 

Of course, the moving targets are affected by time 
domain truncation just as the stationary targets are. Their 
spectra are also spread throughout the passband and stopband 
with the attendant loss of efficiency. 

We have calculated several examples to help provide an 
insight into these truncation effects. The assumptions for 
our calculations are: 

1. Speed of moving target is such that the response 
of the MTI processor is at a maximum. 

2. PRF is fixed at f^. = 1/T. 

3. Pulse width is identical for moving and 
stationary targets and is given by t. 
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Rcjure 3.11. Freqiioncy spectrum of stationary target where echo pulse train consi 
of k = 2 hits per scan 
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Figure 3.12. Frequency spectrun: of stationary target where echo pulse train 
consists of k = 6 nits per scan 
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Fipure 3.13. I’rcquoncy spectrum of stationary target v/hore echo pulse train 
consists of k = 10 hits per scan 



4 . There is equal input energy from stationary 
and moving targets . 

First we calculate the output spectriim based on the 
nonideal input spectrum and the frequency transfer function 
of the filter: 



Y^{f) = H{f) X^{f) (3-9) 

where is the output spectrum for an input X^^ is 

defined as the nonideal input spectrum of a truncated 
stationary target echo train. H(f) is the transfer function 
of the filter. 

Also we have, for moving targets: 



Y2(f) = H(f) X2(f) 



(3-10) 



Following this calculation, we form the integral 






1/T 



|Y^(f) p df 



(3-11) 



and 






1/T 



(3-12) 



|Y2(f) 1 df 



Finally, we form a ratio of the above quantities. Let 



R = 



E ' 
^1 



^db = ^ 



(3-13) 
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Thus R is the measure of energy contained in the output from 
an input, truncated, moving target echo train relative to 
the energy contained in the output caused by a truncated, 
stationary target echo train, with the four assumptions 
made above . 

In the ideal case, this ratio would be infinite. In 
the nonideal case (nonideal, yet still tractable since the 
behavior can be predicted) , the ratio will never be infinite. 
The ratio will form what is an upper bound for the radar 
system under these conditions. Under no circumstance, can 
the user expect to see a greater difference between output 
moving target echo energy and output stationary target 
echo energy than R (again, under the aforementioned assump- 
tions) . 

The results of these calculations are tabulated and 
shown in Table 3.1. 

As might be expected from an intuitive approach, the 
larger feedback enables the filter to achieve a higher 
signal-to-noise ratio (where noise is defined as the sta- 
tionary target echo) . Also, as expected, the ratio rises 
with the number of hits on target. There are no real sur- 
prises "here, but it is emphasized, that prior to the design 
of a radar MTI processor, one must decide on many specifica- 
tions. One of the important specifications is the expected 
number of "hits" on target. The problem is tractable either 
for a "least number of hits" expected, or some mean number 
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TABLE 3-1 



ENERGY RATIO : Energy contained in the output of 

an MTI processor caused by a moving target (truncated 
pulse train) compared to the energy contained in the output 
caused by a stationary target (truncated pulse train) . 

This figure represents the maximum ”Signal-to-Noise" ratio 
possible given the assumptions listed in the text. Ratio 
is in decibels. 
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First Order Feedback Coefficient 
0.3 0.7 0.9 



2 6.7 12.3 15.9 



10 15.3 19.3 23.5 



20 18.5 22.9 27.0 
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of hits. But, most importantly, the problem must be faced 
at the beginning of design, not, as in some cases related 
to us, after construction of the radar when the best answer 
is a retrofit which never really works. 

As an aside, it has occurred to us that since we do 
know the "noise” spectrum and signal spectrtam, that it is 
possible to design an optimxim filter (maximizes the signal- 
to-noise ratio) . Such an optimum filter can be designed 
using the derivation given by Taub and Schilling [31] for 
nongaussian noise; 



H^p^(f) = X*(f)/G(f) (3-14) 

where X*(f) is the complex conjugate of the input signal 
spectrum, G(f) is the noise energy spectrum, and H^p^(f) 
is the transfer function for the optimum filter. Figure 
3.14 shows the computer-aided optimum transfer function. 

This transfer function is admittedly very limited. It 
is only optimum for the four assumptions given above. If 
the target changes radial speed, then the filter is no 
longer optimum, a fact which probably precludes its use. 

On the other hand, optimum filters could be designed for a 
whole bank of doppler frequencies (radial speeds) and the 
output of each compared after each input echo. Such a 
scheme, while theoretically sound, is probably more than 
is really needed in the field of MTI . 
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Figure 3.14. Transfer function of cptiiaum filter for moving target whose doppler free 
half the PRF vs. stationary target. 



Section 3.3 The Recursive Comb Filter as the MTI 
Processing Element 

The original AN/SPS-58 air search radar utilized a 
digital recursive filter as its main MTI processor. The 
OPEVAL of this radar was conducted in 1972. In 1975, the 
radar set was still undergoing redesign. One of the prob- 
lems associated with this radar was the difficulty that was 
experienced with second-time-around radar targets . A second- 
time-around target is a target (moving or stationary) whose 
echo appears in the pulse repetition interval succeeding 
the one in which it should appear. That is, the echo appears 
after transmission of the succeeding pulse. The echo appears 
as if it were at a small range from the radar, when in 
actuality, it is at a great range. The actual ambiguity 
is equal to the maximum range of the radar. If, for in- 
stance, the echo appears such that the apparent range is, 
say 50 miles, and the maximum range of the radar is 250 
miles, then the actual range is 300 miles. 

One means to eliminate, or at least categorize, second- 
time-around echoes is to change the PRF and see if the echo 
changes its position drastically. Changing the PRF naturally 
changes the maximum range of the radar. If the echo is a 
first-time-around echo — a legitimate echo — then changing 
the PRF will not affect its position on the PPI scope. How- 
ever, if it is a second-time-around echo, then changing the 
PRF such that the maximum range of the radar is, say 270 
miles, will shift the position of the echo by 20 miles so 
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that its apparent range is 30 miles. When the target "jximps" 
like this, the operator is trained to conclude that his 
target is really a second-time-around echo. 

The OPEVAL of the AN/SPS-58 (Westinghouse, prime con- 
tractor) continued to have severe problems with second-time- 
around echoes into 1975. The reason given for NOT providing 
the operator with a PRF dither was 

"Strong land clutter does break through at ' 
the top of the display at a range of about 
19 miles but the MTI performance is very 
good. The recursive digital MTI is not 
amenable to quick PRF changes so that PRF 
dither and range correlation for the 
removal of second-time-around targets is 
not a feature of this radar."* 

The reason for believing that one could not change the 
PRF with a recursive digital filter is a long-standing 
tradition among radar designers, who believe that the 
length of the transient is sufficient to interfere with 
cancellation, and that it is not desirable to provide PRF 
changes in a recursive filter. Part of this feeling stems 
from a reading of Skolnik's excellent text. In it, he warns 
that second-time-around echoes cannot be cancelled by PRF 
changes ... that, as a matter of fact, staggered PRF's 
should be avoided if second-time-around echoes from sta- 
tionary targets are a problem. His analysis is thorough and 



* 

Private communication with Mr. James Whitaker, NELC, 
Code 1260 who quoted official documents dealing with the 
AN/SPS-58 OPEVAL of 18 July 1972. 
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correct. However, the interpretation on his words is some- 
times inaccurate. For instance, there is NO reason at all 
to avoid changing the PRF with a recursive filter so long 
as the timing circuits in the MTI processor are adjusted so 
that one always compares echoes at the same range. This 
synchronization must be performed anyway (no matter what 
form of filter is used) whenever PRF "dither” or staggering 
is done. This concept can become clearer by the following 
example. Figure 3.15 shows a transmitted pulse train where, 
for simplicity, between each pair of transmitted pulses, 
there are ten range gates . These range gates occupy a 
certain range band (if the maximum range is 100 miles, then 
each range gate contains the energy received from a range 
band of 10 miles) . 

After the second transmitted pulse, the maximum range 
is changed to 80 miles, but ten range gates still are present. 
Now, by simple division, we find that the range band width 
is only 8 miles . 

Since there are ten range gates, we shall be using an 
imaginary filter with ten delay elements (a 10th order fil- 
ter) . After the first pulse is transmitted, echoes (or 
lack thereof) for the ten range gates are received and 
stored successively in the delay elements. After the second 
pulse is transmitted, the contents of these range gates 
are consecutively compared with the new information (the 
echoes from the second pulse) and we get, serially, ten 
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transmitted pulses 




Figure 3.15. Transmitted pulse train with staggered PRF. Ten range 
gates of information introduced to MTI after each pulse. 



transmitted pulses 




range gates 

Figure 3.16. Transmitted pulse train with staggered PRF. Eight range 
gates (of equal range-width) of information introduced to MTI after 
each pulse. The farthest ranges are not processed by MTI. 



297 



outputs ... one from each range gate. In the sixth range 
gate, for instance, we are comparing data from the range 
band of 50 miles to 60 miles. The output is a function of 
the two echoes received from that band of ranges. 

Now the pulse repetition interval is changed for the 
third pulse. There are still ten delay elements and ten 
range gates, but each range band is only eight miles wide. 
The output from the filter as the echoes from the third 
pulse arrive no longer makes any sense. The reason is that 
when we compare the contents of the sixth range gate with 
the arriving echo, we are comparing energy in the 50-60 mile 
range to energy in the 40-48 mile range for the third pulse. 
No conclusions can be drawn from such a comparison. 

But this type of misalignment can exist in transversal 
filters too. Depending on the order of the filter, the 
effect is either serious or tolerable. It is well to be 
reminded that the longer the transversal filter is, the 
more closely it resembles a recursive filter. 

This alignment is rather easily avoided. Most of the 
time, the MTI processing is not required on range bands far 
from the transmitter. In this case, it is a relatively 
easy engineering task to synchronize the filter to begin 
each processing phase after each pulse, and to maintain 
range gate width no matter what PRF is used. Figure 3.16 
shows the conceptual pulse train for this plan. Notice that 
there are ranges for which no processing in the MTI takes 
place. This turns out to be acceptable in most applications. 
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